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odd-order  networks  as  well  as  that  of  several  specific  odd-order  topolo¬ 
gies  of  special  interest  are  discussed.  Practical  monolithic  input*  out¬ 
put*  and  interstage  matching  network  topologies  using  both  lumped  and  dis¬ 
tributed  elements  are  shown,  and  the  potential  for  reducing  required  chip 
area  by  using  networks  of  3rd  order  in  lieu  of  4*^  order  is  shown  for  both 
low-noise  and  power  amplifier  examples. 
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Chapter  1 


Introduction 


While  microwave  integrated  circuits  of  hybrid  construction  allow  con¬ 
siderable  modifications  to  a  design  to  optimize  performance  after  fabrica¬ 
tion.  monolithic  microwave  integrated  circuits  (MMIC's)  make  this  impossi¬ 
ble  in  most  cases.  For  this  reason,  accurate  design  and  fabrication  of 
MMIC's  is  of  much  greater  importance  than  ever  before  with  other  fabrica¬ 
tion  technologies. 

Of  central  importance  in  the  MMIC  design  process  is  the  availability 
of  computer  programs  to  assist  the  circuit  designer  in  active  device 
modelling,  synthesis  of  impedance  matching  networks  with  precisely  con¬ 
trolled  responses  for  amplifier  gain  equalization  over  arbitrary 
bandwidtha.  and  amplifier  performance  analysis  and  optimization  by  simula¬ 
tion  prior  to  committing  a  design  for  fabrication.  This  work  extends  the 
state  of  the  art  in  aeveral  of  these  areas,  and  includes  original  computer 
programs  that  simplify  and  make  more  convenient  some  of  the  the  tasks  that 
must  be  accomplished  during  the  design  of  monolithic  microwave  amplifiers. 

Chapter  2  extends  the  exact  lumped  element  matching  network  synthesis 
procedures  of  Petersen  [1]  (valid  only  for  even-order  gain  functions)  to 
allow  exact  synthesis  of  equiripple  lumped  element  networks  of  odd  order 
as  well.  The  unique  traits  of  these  odd-order  networks  are  examined,  and 
procedures  for  taking  maximum  advantage  of  these  during  the  synthesis  pro¬ 
cess  are  presented.  Examples  are  used  to  illustrate  the  significant 
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advantages  in  circuit  simplicity  and  reduced  physical  size  now  available 
to  the  designer  through  the  use  of  3rd  order  networks  in  many  applications 
where  4**1  order  matching  networks  were  previously  required. 

Chapter  3  extends  the  theory  for  exact  synthesis  of  distributed  ele¬ 
ment  matching  networks  to  include  odd-order  networks.  The  unique  traits 
of  these  distributed  element  odd-order  networks  are  also  examined.  List¬ 
ings  of  the  computer  subroutines  used  to  implement  the  new  gain  function 
approximation  algorithms  for  both  the  distributed  and  lumped  element  cases 
are  presented  in  Appendix  A. 

The  derivation  of  an  analytical  method  for  determining  the  theoreti¬ 
cal  limits  of  the  new  synthesis  parameter  available  with  odd-order  func¬ 
tions  (the  denominator  polynomial  constant)  is  presented  for  3rd  order 
networks  in  Chapter  4.  The  procedure  is  completely  general  in  that  it 
applies  to  all  3rd  order  networks  —  regardless  of  bandwidth*  ripple* 
minimum  insertion  loss*  or  gain  slope  specifications.  Analytical  methods 
are  presented  for  both  the  lumped  element  and  distributed  element  cases. 
The  computer  implementation  of  these  techniques  is  presented  in  Appendix 
A. 

Chapter  5  presents  design  techniques  for  matching  networks  allowing 
only  cascaded  transmission  lines  and  shunt  open-circuited  stubs  as  circuit 
elements.  These  networks  greatly  simplify  the  fabrication  of  monolithic 
amplifier  circuits*  and  the  advantages  and  limitations  of  this  restricted 
topology  class  of  networks  are  discussed.  Design  examples  for  octave  and 
10Z  band  monolithic  FET  amplifiers  are  presented;  the  measured  RF  perfor¬ 
mance  of  the  octave  band  amplifier  is  also  presented. 
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Tables  of  circuit  element  values  for  several  example  lumped  element 
matching  network  topologiea  are  presented  in  Appendix  B  for  both  1SZ  and 
100Z  bandwidths  as  well  as  several  different  gain  slopes.  Appendix  C 
presents  similar  tables  for  several  practical  distributed  element  network 
topologies. 

Appendix  D  presents  a  concise  summary  of  single-frequency  matching 
techniques,  extremely  simple  lumped  element  circuits  that  may  be  used  to 
implement  these  for  input,  output,  and  interstage  networks,  and  a  computer 
implementation  of  these  procedures  for  convenience. 
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Lap«4  limit  >roi4b»4  litchiig  letvork  Design 
(Arbitrary  Order) 


Computer  programs  have  been  available  for  several  years  for  the  syn¬ 
thesis  of  amplifier  impedance  matching  networks  with  even-order  bandpass 
gain  responses  [1-2].  A  method  for  synthesizing  odd-order  low-pass  net¬ 
works  has  also  been  reported  recently  [3]. 

The  odd-order  synthesis  method  discussed  by  Riddle  and  Trew  in  [3] 
is  based  upon  the  use  of  even-order  low-pass  gain  functions  to  define 
network  response,  and  then  adding  a  multiplicative  factor  involving  the 
complex  variable  to  the  resulting  8jj(p)  function  so  that  an  additional 

circuit  element  can  be  extracted  during  synthesis.  Although  the  authors 
only  discussed  this  technique  in  the  context  of  low-pass  networks,  it 
could  also  be  applied  to  bandpass  gain  functions.  The  odd-order  syn¬ 
thesis  techniques  described  in  this  work  are  fundamentally  different  from 
those  of  Riddle  and  Trew. 

This  chapter  presents  a  technique  for  synthesizing  lumped  element 
bandpass  impedance  matching  networks  with  arbitrary  gain  slope  from  ndd- 
order  gain  functions.  This  is  accomplished  by  extending  the  even-order 
synthesis  procedure  so  that  odd-order  functions  may  be  included  while 
maintaining  realizability  in  the  gain  function.  The  procedure  can  pro¬ 
duce  narrowband  or  wideband  gain  responses  with  equal  facility. 
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A  primary  concern  in  the  design  of  Monolithic  Microwave  Integrated 
Circuits  (MMIC's)  is  to  find  impedance  matching  networks  for  the  active 
devices  that  require  a  minimum  number  of  elements  (for  minimum  wafer 
area)  as  well  as  providing  a  means  of  introducing  FET  bias  voltages  with 
the  fewest  possible  additional  RF  chokes  and  DC  blocking  and  bypass  capa¬ 
citors. 

All  low-pass  matching  networks*  since  they  are  composed  only  of 
series  inductors  and  shunt  capacitors*  require  the  addition  of  a  physi¬ 
cally  large  shunt  inductor*  one  or  more  bypass  capacitors*  and  a  blocking 
capacitor,  before  biasing  provisions  can  be  included  on-chip;  these 
requirements  largely  eliminate  low-pass  networks  from  GaAs  monolithic 
applications.  Bandpass  networks,  however*  allow  shunt  inductors  and 
series  capacitors,  so  it  is  possible  to  employ  bandpass  topologies  that 
only  require  the  addition  of  a  single  DC  bypass  capacitor  to  allow  on- 
chip  biasing.  This  feature  makes  bandpass  networks  the  optimum  choice 
for  MMIC's. 


2..JL.  Synthesis  Algorithm 

Lumped  element  broadband  matching  networks  with  bandpass  gain 
responses  of  arbitrary  gain  slope  can  be  synthesized  by  the  following 
procedure: 

(1)  Specify  a  circuit  topology  (which*  in  turn*  fixes  network  order). 

(2)  Specify  the  desired  bandwidth,  gain  slope*  minimum  insertion 
loss*  and  gain  ripple  of  the  network. 


(3)  Fora  an  equiripple  approximation  [G(x)3  to  the  desired  gain 
response.  (Find  unknown  coefficients  of  G(x)  by  computer  algorithm.) 

A 

(4)  Transform  G(x)  to  G(«  ),  a  transducer  power  gain  function,  by 

2 

substituting  u  *  x. 

(5)  Compute  Jsj^pJp.  then  factor  into  s^p)  *nd  where 

s^tp)  includes  all  the  left  half-plane  poles  and  any  distribution  of  the 
zeroes  of  kn(p)f- 

(Use  js11(p)j2  =  1  -  1*21^^*  wher*  P  =  the  co“Plex  variable.) 

(6)  Compute  z.j(p)  from  Sj^p)  by  using 

1  +  «n(p) 
zll(p>  =  1  -  su(P) 

(7)  Extract  network  element  values  by  continued  fraction  expansion  of 

«n<p>. 

The  following  sections  will  discuss  in  detail  the  solution  to  the 
gain  response  approximation  problem  [step  (3)]»  and  a  text  on  network 
synthesis,  such  as  Weinberg  [4],  may  be  consulted  for  more  detail  on 
steps  (5)  through  (7). 

2.2.  Motivation  for  Odd  Order 

The  techniques  developed  by  Petersen  [1]  can  be  used  to  synthesize 
even-order  bandpass  networks  (i.e..  K=2.4.6...)»  and  these  have  been 
implemented  in  the  CADSYN  commercial  computer  program  [2].  The  simplest 
even-order  network  that  can  be  expected  to  simultaneously  absorb  an  FET's 


parasitic  capacitance  while  providing  a  specific  isipedance  transformation 
is  of  4th  order • 

Figure  2.1(a)  illustrates  a  typical,  practical  4th  order  bandpass 

network  that  can  be  used  to  exactly  absorb  the  parasitic  output  capaci- 

tance  of  an  FET  (C  )  while  providing  an  impedance  transformation  between 
P 

the  FET  output  resistance  (RQut)  and  the  amplifier  load  as  well  as  adding 
gain  slope  compensation  to  provide  flat  gain  across  the  passband  despite 
the  FET's  intrinsic  gain  roll-off  with  frequency.  After  adding  one  capa¬ 
citor.  the  circuit  in  Figure  2.1(b)  can  offer  essentially  the  same  RF 
performance  as  that  in  Figure  2.1(a).  but  with  on-chip  biasing  as  well. 

We  would  like  to  eliminate  two  of  the  inductors  in  the  network  of 
Figure  2.1.  and  use  a  3rd  order  network  such  as  that  shown  in  Figure  2.2 
to  still  provide  simultaneous  exact  parasitic  absorption,  exact  impedance 
transformation,  and  bias  introduction,  while  requiring  the  fewest  possi¬ 
ble  circuit  elements. 

2*1.  Even-Order  Equiripple  Gain  Function  Syathcaia 

The  function  shown  in  equation  (2-1)  will  be  used  to  approximate  a 
specified  ideal  gain  response. 

«*>  •  y*r  <2-» 

where:  x  =  the  frequency  squared  variable  (u^) 

k  =  number  of  high-pass  elements  in  the  network 
(determined  once  topology  has  been  specified) 
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(a) 


Figure  2.2. 


3r<*  Order  Bendpeee  Networks 
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K  =  gain  constant  chosen  to  maintain  realizability 
(0  £  G(x)  £  1  V  x  >  0  ) 

Pg(x)  =  an  Nth  order  polynomial  in  x»  whose  coefficients 
are  to  be  determined* 


An  equiripple  approximation  will  be  computed  by  solving  for  the  unk¬ 
nown  K  and  the  coefficients  of  Pu(x). 

Our  synthesis  procedure  requires  the  initial  specification  of  a 
topology*  as  well  as  bandwidth*  ripple*  gain  slope*  and  minimum  insertion 
loss  for  the  network;  the  algorithm  then  proceeds  to  solve  for  the  coef¬ 
ficients  of  ?N(x)  bo  that  the  specified  approximation  is  realized. 

The  N+l  unknown  coefficients  of  PN(x)  can  be  determined  (for  even  N) 

by  setting  up  a  system  of  N+l  equations  involving  eq  (2-1)  with  the  gain 
value  (G(x) )  specified  at  N+l  critical  points  across  the  passband  alter¬ 
nating  above  and  below  the  specified  ideal  gain  response*  as  shown  in 
Figure  2.3  for  N  =  4.  He  require  G(x)  to  be  decreasing  as  x  increases  or 
decreases  away  from  the  passband  (so  0  S  G(x)  £  1  V  x  >  0);  this  formu¬ 
lation  facilitates  this  behavior  in  the  function  for  even  N»  since  we  can 
specify  an  odd  number  of  critical  points. 


The  ideal  response  shown  in  Figure  2.3  will  be  referred  to  later  as 
the  reference  gain  function*  G^^tx)*  defined  as 


Gref(x) 


10 


r 

I 


i 


GS.  *  _*iq1x)  ] 

10  log10(2)  j 


where  GS  is  the  specified  gain  slope  in  dB  per  octave 


with  passive  elements  only  -  which  is  the  sole  case  of  interest*) 

This  problem  can  be  avoided  by  specifying  the  approximating  function 
with  only  N  critical  points  when  N  is  odd  >  thus  forcing  the  required  even 
symmetry  in  G(x)  and  allowing  s  realisable  gain  function  to  be  generated* 
This  approach  produces  a  system  of  N  simultaneous  linear  equations*  which 
allows  determination  of  N  of  the  N+l  unknown  coefficients  of  P^tx) .  If 
we  choose  one  of  the  coefficients  of  PN(x)  (arbitrarily,  the  constant 

term)  to  be  a  parameter  of  the  synthesis  (i.e.*  chosen  initially,  and  not 
allowed  to  vary) •  then  we  can  solve  the  approximation  problem  with  an 
odd-order  G(x) . 

2*A*1*  Overview  of.  Appro* im«t  ion  Algorithm  -  Odd.  Order 

The  approximation  problem  is  well  formed  once  the  bandwidth,  gain 
slope,  and  minimum  insertion  loss  of  the  ideal  gain  function  as  well  as 
the  allowable  ripple  radius  of  the  approximation  have  been  specified. 

The  computational  task  is  then  to  determine  the  necessary  denominator 
polynomial  coefficients  of  G(x)  to  obtain  the  specified  approximation. 

This  problem  is  solved  by  using  several  iterations  of  a  two-step 
procedure,  once  a  value  for  the  constant  term  of  ?N<x)  has  been  chosen: 

(1)  Assume  locations  for  the  points  of  maximum  error  between 
the  specified  ideal  response  and  the  approximating  function  (the 
critical  points),  and  compute  required  values  of  G(x)  at  these 
points.  Use  this  information  to  specify  a  system  of  N  linear 


u 


(FREQ2) 


Figure  2.4.  Unrealizable  Odd-Order  Gain  Functions 
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equations  in  N  unknown  coefficients  of  P^(x);  solve  the  system  for 
the  N  coefficients. 

(2)  Using  the  PN(x)  coefficients  determined  above*  examine 
G(x)  to  determine  the  true  locations  of  the  critical  points.  Do 
this  by  solving  a  one-dimensional  non-linear  equation  with  multiple 
solutions  via  the  Newton-Raphson  Method  [5].  Using  these  new 
values  for  the  critical  points*  re-enter  step  (1). 

Initially*  our  assumed  critical  points  (used  to  solve  for  the  coef¬ 
ficients)  vill  be  in  error  when  compared  to  the  true  critical  points 
(determined  after  finding  the  coefficients).  If  these  new  critical 
points  are  now  used  to  again  solve  for  the  polynomial  coefficients*  a 
G(x)  will  be  obtained  whose  true  critical  points  are  much  closer  to  the 
values  used  to  compute  the  coefficients.  This  iterative  process  con¬ 
verges  very  rapidly  (usually  within  4  iterations)*  and  is  terminated  when 
the  values  for  the  critical  points  used  in  determining  the  coefficients 
are  found  to  be  within  a  specified  tolerance  of  the  actual  values  for  the 
critical  points  found  when  subsequently  examining  the  resulting  G(x) . 

Besides  incorporating  odd-order  functions*  this  method  differs  from 
Petersen's  Cl]  in  the  method  of  solving  the  non-linear  portion  of  the 
problem  for  the  critical  points.  Where  Petersen  dealt  with  it  as  an  N- 
dimensional  non-linear  problem  requiring  computation  of  a  Jacobian  to 
solve*  this  method  treats  it  as  a  one-dimensional  non-linear  problem  with 
multiple  solutions. 


IS 


2.A>2*  Details  of.  the  Appr^Mmatinn  Algorithm  -  QM  Order 


The  following  procedure  is  used: 


(1)  Baaed  upon  specified  network  order  and  bandwidth*  assume  N  crit¬ 
ical  frequencies  (x^»  i=l»2....*N)  of  alternating  error  (relative  to  the 

ideal  function)  across  the  passband. 

(2)  Compute  required  gain  level*  G^e^(x^).  *t  each  critical  point 
based  on  gain  response  specifications*  assuming  for  the  moment  that  the 
numerator  gain  conatant*  K*  is  unity.  (K  will  be  set  later  to  insure 
that  the  maximum  value  of  G(x)  V  x  >  0  is  1 .  For  the  moment*  no  general¬ 
ity  is  lost  by  assuming  K=1 .) 

(3)  Now  we  have 


G 

req 


(x^ 


Vxi> 


1  2  1  »2  MSS  |N 


(2-2) 


P»<xi)  *  *0  *  J/j  *  xiJ 


aft  ♦  5  a  .  *  x . J 
0  j=l  J  1 


G  (x.)  ’ 
req  i 


i  =  1  *2  * .  •  •  *N 


(2-3) 


The  x^'s  «re  assumed  known  at  this  point*  so  the  N  values  for 

G  (x.)  can  be  computed.  The  coefficients  a.*  j  =  0*1*2 . N  are  to  be 

req  l  r  J 

determined. 


Equation  (2-3)  represents  a  system  of  N  linear  equations  in  N+l  unk¬ 
nowns.  To  allow  solution  of  this  system*  sq  is  chosen  to  be  a  p«i,«— 
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rather  than  a  variable*  and  a  (positive)  value  for  it  is  chosen  at  the 
start  of  the  procedure;  this  reduces  eq  (2-3)  to  N  equations  in  N  unk¬ 
nowns. 

(4)  Equation  (2-3)  may  be  re-vritten  in  the  form 


where : 


X  a  =  b 


N  N-l  2 

tj  Xj  .  Xj  Xj 

N  N-l  2 

l2  x2  *  x2  x2 


(2-4) 


(2-5) 


N  N-l 
*N  XN 


2 

XN  ^ 


=  [  V  *N-1 . *2'  *1  jT 


(2-6) 


I  k  k 

'  f  Greq(xl}  0  ’  Gr«q(x2> 


a0’  *  *  * 


*  Greq(xN>"‘°  ] 


(2-7) 


For  best  computational  speed  and  numerical  stability*  eq  (2-4)  is 
solved  for  a  by  LU  factorization  of  the  X  matrix  [6]. 


factor 

X 

=  LU 

(product  of  triangular  matrices) 

then 

LU  a 

=  b 

from  eq  (2-4) 

solve 

L  7 

=  b 

(a  triangular  system) 

solve 

U  a 

=  c 

(a  triangular  system) 

(In  the  computer  program*  these  steps  are  handled  by  calls  to  subroutines 
DGEFA  and  DGESL  from  the  LINPACK  subroutine  library  [7].) 
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It  would  also  be  possible  to  solve  eq  (2-4)  by  inverting  the  X 
matrix *  and  since  X  is  of  Vandermonde  form*  thia  inverse  could  be  com¬ 
puted  with  less  than  the  normally  required  amount  of  computation*  Since 

2  3 

it  requires  3N  work  to  invert  a  Vandermonde  matrix  and  N  / 3  work  to  com¬ 
pute  the  LU  decomposition*  the  LU  decomposition  is  still  the  most  effi¬ 
cient  way  to  solve  eq  (2-4)  in  practical  problems*  however*  where  6**1 
order  is  the  largest  network  usually  considered. 


[  NOTE;  At  this  point*  the  N+l  denominator  polynomial  coefficients  of  eq 
(2-3)  have  been  determined.  The  next  major  task  is  to  improve  bur  esti¬ 
mate  of  the  locations  of  the  x^,  by  examining  the  current  G(x).  and  find¬ 
ing  its  true  critical  points.  The  remainder  of  the  algorithm  deals  with 
finding  these  critical  points.] 


(5)  If  Gref(x)  the  specified  ideal  (reference)  gain  function  that 

is  being  approximated*  then  the  critical  points  internal  to  the  passband 
are  those  points  (x^)  where 


d_r 

dxl 


G(x) 


1  =  Jj-T  G  f(x) 
J  dx^  ref 


1 

J 


(2-8) 


The  two  critical  points  at  the  passband  edges  are  fixed*  and  therefore 
known,  throughout  the  procedure.  Ve  need  to  determine  only  the  M-2  crit¬ 
ical  points  internal  to.  the  passband  by  solving  eq  (2-8). 


Ve  expect  to  find  M-2  solutions*  since  we  have  specified  the  approx- 
imating  function  to  have  N-2  critical  points  internal  to  the  passband. 

By  expanding  eq  (2-8)  we  obtain 
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d_r 

dxl 


1 


■*«'  V«> ! 


_  dx 


a 

h*L 


r 


(x) 


-  d-Tp  (X)l 
dxt  M  -  J 


•  m 

*  x 


P„2<,> 


£  Sr.f<*>  ] 


(2-9) 


In  eq  (2-9)  the  coefficient*  of  Pjj(x)  ere  considered  fixed  at  the 

values  determined  for  them  in  step  (4),  so  eq  (2-9)  represents  a  non¬ 
linear  equation  in  one  unknown  (x)»  with  N-2  solutions.  One  good  way  to 
find  these  solutions  is  via  the  Mewton-Raphson  Method. 

(6)  We  wish  to  find  all  x  across  the  passband  such  that 


-  £[«,.£<*>]  ■  0  * 

Equation  (2-10)  defines  the  objective  function.  F(x),  which  we  will  seek 
to  minimize  when  finding  the  critical  frequencies.  We  know  that  the 
first  and  second  derivatives  of  G(x)  exist  because  it  is  a  rational  func¬ 
tion  of  polynomials;  the  derivatives  of  Gref(x)  exist  because  it  is  an 

exponential  function.  Then,  starting  from  an  initial  estimate  for  one  of 
the  solutions.  x.^k^  (the  ith  critical  point),  we  can  compute  an  impriw»H 

estimate  for  that  solution.  x.^k+*\  via  the  Newton-Raphson  Method  and 


(k+1) 


F(x.(k)) 

<LJF(x  00)1 
dxlFUi  'j 


p  S  1  (2-11) 


Where  p  is  a  parameter  used  to  reduce  the  step  size  in  the  early  phases 
of  the  Newton  algorithm.  The  update  specified  in  eq  (2-11)  is  iterated 
until  some  stopping  criteria  is  satisfied. 
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One  property  of  this  method  is  its  excellent  convergence  rate  once 

(k)  # 

x.  approaches  "close  enough"  to  the  true  solution*  x^  ;  another  is 

(k)  * 

its  quite  unpredictable  performance  when  x^  is  "too  far"  from  x^  • 

This  second  property  makes  it  necessary  to  take  special  precautions  in 
using  this  method  in  an  algorithm  to  find  all  the  solutions  to  eq  (2-9)* 
(7)  During  the  first  few  overall  iterations  of  the  approximation 
procedure*  the  initial  estimates  of  the  locations  of  the  critical  points 
are  typically  not  very  good.  This  means  that  these  estimates  cannot  be 
used  as  reliable  starting  points  for  the  Newton  iterations  to  find  the 
critical  points  for  two  reasons:  (1)  the  Newton  Method  would  sometimes 
converge  to  a  point  outside  the  desired  passband*  and  (2)  the  method 
might  sometimes  "miss"  a  critical  point  by  converging  to  a  single  point 
from  two  different  starting  points,  thus  producing  erroneous  results  (all 
critical  frequencies  will  not  be  found) .  An  alternate  strategy  was  dev* 
ised  to  avoid  these  problems. 

During  the  first  few  overall  iterations*  the  Newton  starting  point 
for  each  of  the  N-2  critical  points  is  chosen  by  a  short  search  prior  to 
entering  the  Newton  iteration.  During  this  search  the  second  derivative 
of  F(x)  is  evaluated*  and  a  point  is  accepted  as  a  starting  point  once 
the  second  derivative  at  that  point  is  non-zero  and  of  the  proper  sign 
for  the  region  of  the  curve  (  G(x)  )  where  the  critical  point  is  being 
sought.  (This  insures  that  the  Newton  iterations  will  proceed  in  the 
proper  direction  from  a  chosen  starting  point*  and  that  we  never  try  to 
evaluate  eq  (2-11)  when  the  denominator  might  evaluate  to  zero.) 

Depending  upon  the  magnitude  of  the  second  derivative*  and  the  dis¬ 
tance  away  from  x*»  the  step  size  computed  by  eq  (2-11)  may  be  much  too 
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large,  to  control  this  problem,  the  update  expression  (based  on  eq 
(2-11))  is  modified  to  allow  a  maximum  step  size  of  only  0*03  at  any 
time.  This  eliminates  erratic  behavior  of  the  algorithm  when  far  from 
the  true  solution,  while  not  hampering  rapid  convergence  once  the  vicin¬ 
ity  of  the  solution  has  been  reached. 

After  the  first  few  iterations  of  the  overall  procedure,  the  loca¬ 
tions  of  the  critical  points  change  very  little  between  iterations,  and 
the  critical  points  found  during  the  previous  iteration  are  used  as  the 
starting  points  for  the  next.  This  minimizes  the  number  of  Newton  itera¬ 
tions  required  once  the  danger  of  erratic  behavior  is  past* 

(8)  After  the  true  critical  points  have  been  determined,  the  overall 
procedure  is  repeated  by  re-entering  at  step  (4).  Once  the  values  for 
the  critical  points  have  been  determined  well  enough  that  aa.  Newton 
iterations  are  required  to  find  any  critical  point  after  re-computing  the 
coefficients,  the  loop  is  exited  with  the  final  values  for  the  coeffi¬ 
cients  and  the  critical  frequencies. 


2.4.1*  Implementation 

The  above  algorithm*  with  very  slight  modification,  can  solve  the 
even-order  equiripple  gain  function  synthesis  problem  as  well.  This 
method  has  been  implemented  in  a  computer  subroutine  for  solving  the 
lumped  element  approximation  problem  with  functions  of  arbitrary  order, 
and  combined  with  other  existing  subroutines  to  form  a  complete  lumped 
element  network  synthesis  computer  program  (designated  ARBSYN).  A 
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commented  FORTRAN  listing  of  the  approximation  subroutine  (designated 
ODDSYN)  is  presented  in  Appendix  A. 


Z.JL*  Effecta  of.  Varying  the  D^nmiminnr  Polynomial  Constant  (DPC) 

While  the  denominator  polynomial  constant  (DFC)  may  be  chosen  as  a 
synthesis  parameter  when  designing  odd-order  networks.  there  is  always  a 
finite  allowable  range  for  CPC  that  will  result  in  a  realizable  gain 
function. 

The  fundamental  realizability  criteria  for  G(x)  is  0  £  G(x)  £  1 
V  x  >  0.  We  can  always  insure  that  max[  G(x)  ]  s  1  by  properly  choosing 
the  numerator  gain  constant*  K,  in  eq  (2-1).  But  the  denominator  polyno¬ 
mial  coefficients  control  the  0  £  min[  G(x)  ]  portion  of  the  criteria. 

To  have  0  £  G(x)  for  very  small  values  of  x*  as  x-*0*  we  need  the 
constant  term  of  PN(x)  to  be  positive.  To  have  0  £  G(x)  for  very  large 

values  of  x.  as  x-*oo.  we  need  the  high  order  coefficient  of  PN(x)  to  be 
positive  as  well. 

In  general*  then*  we  can  say  that  the  (open)  interval  of  allowable 

it  ^ 

DPC  values  is  (  O.DPC  )•  where  DPC  is  the  smallest  positive  value  that 
produces  a  P^(x)  with  a  negative  high  order  coefficient.  While  the  value 
of  DPC*  varies  widely  with  the  gain  function  specifications  being  approx¬ 
imated*  it  is  usually  on  the  order  of  0.01  to  1.0  for  typical  networks. 

Any  DPC  value  chosen  from  the  realizable  interval  (  O.DPC*  )  will 
result  in  a  gain  function  that  meets  the  approximation  specifications  and 
allows  synthesis  of  a  physical  network  corresponding  to  that  gain 


function*  but  the  circuit  element  values  and  required  terminating 
impedance  can  vary  quite  markedly  as  DPC  is  al loved  to  range  within  this 
realizable  interval. 


2.1.1.  Gain-Bandwidth  Behavior  with  Q££ 

One  of  the  more  important  performance  statistics  of  a  matching  net¬ 
work  is  its  gain-bandwidth  product  (GBW)*  and  for  odd-order  networks  GBW 
can  be  observed  to  vary  considerably  with  DPC.  The  normalized  gain- 
bandwidth  product  of  a  network  is  directly  proportional  to  the  value  of 
(the  first  capacitive  element)  when  Cj  is  connected  in  shunt*  and  pro¬ 
portional  to  1/Cj  when  is  connected  in  series.  Using  this*  we  can 
examine  the  GBW  of  several  different  networks  as  DPC  is  varied  over  its 
realizable  limits. 

Each  of  the  four  3rd  order  networks  shown  in  Figures  2.5  and  2.6 
were  synthesized  for  0  gain  slope*  0.05  dB  ripple*  0  minimum  insertion 
loss*  and  octave  bandwidth.  In  each  case  DPC  was  varied  over  the  full 
realizable  range*  and  the  plots  show  the  resulting  behavior  of  GBW  for 
each  case.  The  lower  limit  on  DPC  is  zero  in  all  cases*  and  the  maximum 
realizable  value  is  determined  by  the  procedures  described  in  Chapter  A. 

In  many  cases*  a  DPC  value  can  be  found  internal  to  the  realizable 
interval  which  maximizes  GBW.  Two  exsmples  of  this  behavior  are  shown  in 
Figure  2.5.  The  value  of  DPC  which  maximizes  GBW  for  one  network*  how¬ 
ever*  is  not  related  to  the  DPC  which  maximizes  GBW  for  another  network* 
nor  are  the  general  "shapes"  of  the  GBW  vs  DPC  curves  necessarily  similar 
for  different  networks. 
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interval  for  DFC.  Figure  2.6  shows  two  examples  where  this  type  of  GBW 
vs  DPC  behavior  is  observed.  In  Figure  2.6(a)*  GBW  can  always  be 
increased  by  reducing  DFC*  though  below  about  0.001  the  effect  is  much 
less  dramatic.  In  Figure  2.6(b)*  with  a  different  network*  GBW  can  be 


Figure  2.6.  Gain-Bandwidth  Effects  of  DPC  Variation 
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increased  by  increasing  DPC  until  the  limiting  value  for  realizability  is 
encountered  (at  about  0.023  in  this  example).  (Mote  that  the  DPC  axis  in 
Figure  2.6(b)  is  not  drawn  to  scale.) 

2.1.2.  Terminating  impedance  Behavior  with  Q££ 

Different  values  for  DPC  in  an  odd-order  network  produce  different 
terminating  resistance  values  in  the  synthesized  network.  Figure  2.7 
shows  two  typical  examples  of  this  behavior,  where  both  networks  were 
designed  for  0  gain  slope.  0  minimum  insertion  loss,  and  0.05  dB  ripple 
over  an  octave  band. 

The  fact  that  choice  of  DPC  affects  R  .  as  well  as  the  size  of  the 

out 

first  reactive  element  in  an  odd-order  matching  network  is  important  in 
the  design  of  practical  matching  networks,  where  a  simple  3rd  order  net¬ 
work  can  often  be  designed  to  absorb  an  FET's  parasitic  capacitance  and 
provide  a  required  impedance  transformation  instead  of  using  a  larger  4^ 
order  network. 

2.1*  Choosing 

The  denominator  polynomial  constant  is  an  additional  synthesis 
parameter,  available  only  when  designing  odd-order  networks,  and  some 
discussion  is  warranted  on  how  best  to  use  it  in  designing  circuits. 


Figure  2.7.  R  .  Effects  of  DFC  Variation 
out 

2.4.1.  simultaneous  Paraaitic  Absorption  and  lapedancr  Matching 

The  aost  common  problem  in  designing  practical  noaoiithic  amplifiers 
is  that  of  finding  Batching  netvorks  with  a  minimum  number  of  elements 
that  allow  simultaneous  absorption  of  sn  FBI's  parasitic  capacitance  and 
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impedance  matching  to  a  specific  load  or  source.  The  best  strategy  for 
this  situation  is  to  simply  use  DPC  in  the  same  manner  that  the  other 
synthesis  parameters  are  used  to  search  out  a  combination  of  slope*  rip¬ 
ple*  bandwidth*  minimum  insertion  loss*  and  DPC  that  allows  solving  the 
matching  problem*  preferably  without  requiring  a  Morton  Transformation 
(see  8  2.12). 

In  this  instance*  the  availability  of  DPC  as  an  extra  synthesis 
parameter  increases  the  likelihood  that  a  suitable  network  of  lower  order 
can  be  found  to  solve  the  particular  matching  problem. 

2.£.2.  Unique  DPC  lot  Maximum  Gain-Bandwidth  Product 

As  discussed  above,  there  is  always  a  DPC  value  (either  internal  to 
the  realizable  interval  for  DPC,  or  at  one  of  the  endpoints)  that  maxim¬ 
izes  the  gain-bandwidth  product  of  a  matching  network.  While  the  optimum 
DPC  for  best  GBW  can  be  found  by  trial  and  error*  a  systematic  search 
routine  has  been  incorporated  into  the  ARBSYN  computer  program  to  allow 
quick  and  convenient  identification  of  the  maximum  gain-bandwidth  solu¬ 
tion  when  synthesizing  odd-order  networks.  This  routine  represents  a 
compromise  between  computer  execution  time  and  the  accuracy  with  which 
the  optimum  DPC  is  identified.  As  implemented,  the  routine  identifies 
the  best  value  for  DPC  within  5Z>  and  does  this  in  30  iterations  or  leas. 

The  search  algorithm  is  as  follows: 

(1)  Initially  examine  the  resulting  GBW  of  networks  synthesized  with 
DPC  values  spanning  the  interval  [0.001*  10.0].  Perform  up  to  10  evalua¬ 
tions.  with  DPC  values  spanning  the  interval  in  a  geometric  progression. 
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Save  the  DPC  value  from  this  sub-search  (DPCj)  that  produces  the  best 
GBW. 


(2)  Seduce  the  search  increment*  and  search  a  much  reduced  interval 
including  DPCj  for  the  value  that  produces  maximum  GBV  (DPC^).  Perform 
up  to  10  evaluations. 

(3)  Again  reduce  the  search  increment*  and  search  the  reduced  inter¬ 
val  which  includes  DPC^  until  the  value  that  produces  maximum  GBW  (DPC^) 
is  found.  Perform  up  to  10  evaluations.  DPC^  is  then  within  5Z  of  the 
true  value  which  optimizes  gain-bandwidth  product . 


Within  each  sub-search*  evaluations  are  performed  by  successively 
increasing  DPC.  Once  a  DPC  value  is  reached  that  produces  an  unrealiz¬ 
able  gain  function*  that  sub-search  is  terminated  and  the  algorithm 
advances  to  the  next  stage.  (Once  an  unrealizable  DPC  has  been  encoun¬ 
tered*  all  larger  values  will  likewise  be  unrealizable.) 

Each  GBW  evaluation  within  a  search  involves  synthesizing  a  network 
with  the  current  DPC  value*  and  noting  the  resulting  GBW  value  from  the 
first  (capacitive)  element  of  the  network. 

Two  other  strategies  were  also  tested  for  selecting  a  unique  DPC* 
neither  of  which  proved  to  be  successful.  The  first  was  to  choose  DPC 
such  that  the  second  derivative  of  the  gain  function  at  the  center  criti¬ 
cal  point  was  zero;  it  turned  out  that  equiripple  solutions  with  zero 
second  derivatives  at  the  center  critical  point  do  not  exist.  The  other 
was  to  choose  DPC  such  that  the  slope  of  G(x)  at  the  upper  passband  edge 
was  made  most  negative*  in  hopes  of  maximizing  gain-bandwidth  perfor¬ 
mance;  it  turned  out  that  this  slope  was  not  correlated  to  gain-bandwidth 
performance . 
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2.1.  Practical  Monolithic  Input  networks 

There  are  only  a  few  input  network  topologies  which  need  to  be  con¬ 
sidered  when  designing  practical  monolithic  amplifiers.  Since  FET  bias 
will  need  to  be  introduced,  we  require  that  the  network  contain  a  shunt 
inductor  as  the  first  element  at  the  FET's  terminals,  so  that  bias  may  be 
introduced  by  simply  adding  a  DC  bypass  capacitor  at  the  ground  end  of 
the  inductor.  We  would  also  like  the  network  to  contain  a  series  capaci¬ 
tor  so  that  no  additional  DC  blocking  capacitor  is  required  at  the  input, 
but  this  is  not  a  necessity.  We  also  require  that  the  4th  order  networks 
contain  an  adjacent  pair  of  series  and  shunt  inductors  so  that  a  Norton 
Transformation  may  be  used,  if  needed  to  match  impedances. 

One  of  the  three  3rd  order  networks  shown  in  Figure  2.8(a-c)  will 

often  be  able  to  handle  any  practical  narrowband  or  wideband  FET  matching 
problem  without  requiring  a  Norton  Transformation,  but  one  of  the  two  4th 
order  networks  at  Figure  2.8(d-e)  will  almost  always  be  able  to  solve  the 
problem  when  the  3rd  order  networks  cannot.  The  4^  order  networks  will 
usually  require  a  Norton  Transformation  to  match  impedances,  and  this  can 
be  accommodated  by  splitting  and  re-arranging  the  two  inductors.  In  the 
figure,  the  assumed  resistive  source  impedance  for  the  amplifier  is  shown 
at  the  left  in  the  networks,  and  the  series  R-C  model  for  the  input  of 
the  FET  is  shown  at  the  right. 

Second  order  networks  are  not  practically  useful  in  themselves, 
since  they  will  almost  never  be  able  to  simultaneously  absorb  an  FET's 
parasitic  capacitance  while  providing  a  required  impedance  transforma¬ 
tion.  Fifth,  sixth,  and  higher  order  networks  are  likewise  not  treated 
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here  (for  input  and  output  networks)  since  they  contain  too  many  circuit 
elements  to  be  of  practical  use  in  MMIC's.  especially  GaAs  MMZC's  where 
the  cost  of  wafer  area  is  so  high. 


2-2..  Practical  Monolithic  Output  Network* 

All  of  the  practical  considerations  discussed  above  apply  equally 
well  to  the  choice  of  a  monolithic  output  network  for  an  FET»  but  the 
topologies  are  slightly  different.  The  circuits  shown  in  Figure  2.9(a-b) 
are  3r<*  order  networks,  while  those  at  Figure  2.9(c-d)  are  order  cir¬ 
cuits.  In  the  figure,  the  shunt  R-C  model  for  the  output  of  the  FET  is 
shown  at  the  left,  while  the  assumed  resistive  load  impedance  for  the 
amplifier  is  shown  at  the  right. 

2.1.  Practical  Monolithic  Interstage  Network* 

The  task  of  designing  a  monolithic  interstage  network  is  almost 
always  much  more  difficult  than  that  of  designing  an  input  or  output  net¬ 
work.  This  is  for  several  reasons:  (1)  the  problem  now  requires  exact 
parasitic  absorption  at  both  ends  of  the  network.  (2)  the  required 
impedance  transformation  (from  the  output  resistance  of  one  FET  to  the 
input  resistance  of  another)  is  usually  quite  extreme,  and  (3)  large 
values  of  gain  slope  compensation  are  often  required. 

We  would  again  like  to  be  able  to  introduce  FET  bias  voltages  with 
the  addition  of  a  minimal  number  of  elements,  so  this  dictates  shunt 
inductors  at  the  FET's  and  a  series  capacitor  to  provide  DC  blocking 


33 


between  the  two  stages  if  possible.  Candidate  topologies  that  meet  these 
criteria  are  shown  in  Figure  2.10.  (The  output  models  for  the  first 
FET's  are  drawn  at  the  left  in  the  figure,  and  the  input  models  for  the 
second  FET's  are  shown  at  the  right.) 

The  4**1  order  interstage  network  at  (a)  in  Figure  2.10  represents 
the  smallest  practical  circuit  that  might  be  synthesized  to  serve  as  an 
interstage  network,  and  it  also  has  the  advantage  that  a  Norton  Transfor¬ 
mation  is  possible  to  assist  in  impedance  matching.  This  network,  how¬ 
ever.  would  require  a  DC  blocking  capacitor  as  well  as  an  additional  RF 
choke  (besides  bypass  capacitors)  to  by  fully  usable. 

The  6***  order  interstage  topology  at  (b)  of  Figure  2.10  is  also  a 
good  candidate  that  requires  only  the  addition  of  bypass  capacitors  to  be 
fully  usable,  and  it  can  potentially  use  a  Norton  Transformation  for 
impedance  matching. 

When  neither  of  the  above  networks  prove  to  be  satisfactory,  there 
are  still  two  more  possibilities.  Neither  can  produce  a  gain  response 
whose  performance  specifications  can  be  predicted  exactly  in  advance 
(unlike  the  others),  but  they  do  present  practical  alternative  solutions 
to  the  interstage  matching  problem. 

In  both  cases  the  general  strategy  is  to  first  form  a  2nd  order  net¬ 
work  at  the  input  to  the  second  FET  by  adding  a  shunt  inductor.  This 
network  should  exactly  absorb  the  FET's  parasitic  input  capacitance,  and 
will  produce  some  (uncontrolled)  impedance  transformation  to  a  resistive 
value  R  .  The  interstage  problem  has  now  been  reduced  to  one  of  matching 
the  FET's  output  to  the  resistive  load  R*—  an  output  matching  problem. 


This  part  of  the  problem  can  either  be  aolved  with  a  3r<*  order  network 
(Figure  2.10(c))  or  a  4th  order  network  (Figure  2.10(d))*  and  the  two 
networks  caacaded  to  form  the  overall  interatage  network. 

Due  to  the  gain  ripple*  and  perhaps  gain  slope*  in  the  responses  of 
the  two  networks,  the  overall  gain  response  of  the  cascaded  network  over 
a  wide  band  is  not  directly  predictable  from  knowledge  of  the  individual 
gain  responses  themselves.  In  practice,  however*  this  is  of  little 
consequence*  since  computer  optimisation  can  be  applied  to  the  cascaded 
network  to  usually  obtain  a  very  satisfactory  gain  response  even  over  an 
octave  band  or  more. 

In  theory*  the  cascaded  technique  is  inferior  to  the  4^.  5^.  or 
6th  order  directly  synthesized  interstage  networks  because:  (1)  it  is 
not  possible  to  exactly  predict  the  overall  gain  response  of  the  cascade* 
(2)  the  cascade  often  requires  some  computer  optimization  to  obtain  a 
satisfactory  wideband  gain  response*  and  (3)  the  gain  response  of  the 
cascade  will  often  not  be  as  "good"  (low  in  ripple)  as  that  of  a 
corresponding  directly  synthesized  network.  In  practical  terms*  however* 
the  cascade  technique  almost  always  provides  a  solution  to  a  difficult 
matching  problem*  even  when  a  directly  synthesized  network  cannot.  It 
can  therefore  be  very  useful  in  designing  MMIC's. 

Regardless  of  the  topology  chosen*  it  is  important  that  all  possible 
distributions  for  the  zeroes  of  the  (p)  function  be  examined  during 

the  search  for  a  suitable  interstage  network.  The  default  distribution 
(all  zeroes  in  the  left  half  complex  plane)  is  often  not  the  optimum  for 
solving  interstage  matching  problems. 


When  using  topologies  (a)  and  (b)  note  that  during  the  initial 


phases  o£  the  synthesis  the  search  is  for  a  network  that  absorbs  Cj 
exactly  with  a  source  resistance  of  Rj  >  and  has  an  RjCj  product  the  same 
as  that  for  the  FET  input  model  being  matched.  After  this  condition  has 
been  achieved  in  a  network,  a  Norton  Transformation  may  be  performed  to 
adjust  the  output  resistance  value  to  that  required  by  the  FET  model. 
Since  the  RC  product  is  invariant  under  impedance  scaling,  the  parasitic 
C  value  will  also  be  correct  once  the  output  resistance  has  been 
adjusted. 


2.1Q..  General  Odd-Order  Gain-Bandwidrh  RefiullA 

It  is  interesting  to  investigate  the  general  gain-bandwidth  behavior 
of  odd-order  networks  and  relate  this  to  the  known  characteristics  of  the 
gain  functions. 

The  fundamental  integral  restrictions  for  gain-bandwidth  originally 
derived  by  Bode  Cl 3]  and  expanded  by  Ku  and  Petersen  [14]  are  useful  for 
examining  how  gain-bandwidth  behavior  is  related  to  the  gain  function. 

For  the  shunt  RC  case  (low-pass  (LP)  constraint)  this  restriction  is 

[  rrk>  ]  s  k  <2-u: 

and  for  the  series  RC  case  (high-pass  (HP)  constraint)  it  is 
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When  comparing  two  networks  for  gain-bandwidth  (GBW)  performance* 
assuming  both  meet  the  same  gain  envelope  specifications  within  the 
passband*  the  "better"  network  is  the  one  with  the  smaller  value  for  the 
GBW  integral  (the  left-hand  side  of  eq  (2-11)  or  (2-12)*  for  the  LP  or  HP 
cases*  respectively).  This  generally  corresponds  to  "steeper  skirts"  in 
the  gain  response  outside  the  specified  passband. 

Figure  2.11(a)  illustrates  why  6th  order  networks  exhibit  better 
gain-bandwidth  performance  than  order  networks  (with  both  designed  to 
the  same  passband  specifications).  The  6^  order  response*  with  3 

passband  ripples*  has  a  sharper  drop  in  gain  both  above  And.  below  the 

passband  than  does  the  4th  order  response  with  2  passband  ripples.  This 

/ 

is  the  familiar  case  for  even-order  bandpass  networks  where  increasing 
network  order  always  produces  a  higher  GBW  product. 

When  synthesizing  a  5th  order  network  to  the  same  passband  specifica¬ 
tions.  however,  something  slightly  different  occurs.  Figure  2.11(b) 
illustrates  the  relative  position  of  typical  4t*'  and  5^  order  gain 
responses.  Since  we  are  still  required  to  synthesize  the  5**1  order  net¬ 
work  with  just  2  passband  ripples  for  realizability*  the  typical  result  is 
a  gain  response  superior  to  the  4th  order  case  at  frequencies  above  the 
passband.  but  worse  than  the  4C^  order  response  below  the  passband. 

To  see  the  reason  for  this*  let's  look  at  what  is  happening  with  the 
denominator  of  the  gain  function  (the  numerators  are  the  same  in  both 
cases).  Figure  2.12  shows  an  approximate  plot  of  typical  4th  and  5th 
order  denominator  polynomial  values  for  a  gain  function  of  the  type  shown 
in  Figure  2.11(b).  At  this  time  it  is  also  helpful  to  recall  the  form  of 
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Figure  2.12.  Typical  Gain  Function  Denominator  Values 

tli  2 

the  5  order  gain  function  (x=u  and  k=number  of  high  pass  elements) 

k 

G(x)  =  T  J — j  2 

OjX  ♦  O^X  ♦  O^X  +  O^X  +  OjX  ♦  Oq 

Above  the  passband  (  x  >  1  )  the  5th  order  denominator  term  causes  the 

denominator  to  grow  positively  at  a  greater  rate  than  that  achieved  by  the 
4^  order  function  —  this  is  expected.  Below  the  passband  (  x  <  Xj  )* 

however*  the  5ttl  order  denominator  cannot  grow  as  rapidly  as  the  4th  order 
function.  This  is  because  since  the  5^  order  function  has  odd  symmetry* 


it  naturally  tends  toward  negative  values  below  the  passband  since  we  syn¬ 
thesise  it  to  tend  toward  increasingly  positive  values  ahove  the  passband. 
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We  can  control  this  tendency  somewhat  by  forcing  the  ordinate  intercept  of 
the  5th  order  function  to  occur  at  some  positive  value  (by  choosing  DPC* 
or  Oq,  to  be  positive)*  but  this  intercept  will  always  be  at  or  below  the 
intercept  of  the  corresponding  4^  order  function.  The  result  of  this 
when  evaluating  the  gain  response  is  that*  below  the  passband*  the  5 
order  gain  response  will  always  be  above  that  of  the  corresponding  4^ 
order  response. 

There  are  two  limiting  cases  where  the  5^  order  performance 
approaches  that  of  the  4**1  order  network.  If  <Xq  is  chosen  so  that  is 
extremely  small,  the  denominator  becomes  effectively  a  4~'r::  order  function* 
and  the  "5th  order"  gain  response  approaches  that  of  the  4th  order  network 
both  above  and  below  the  passband.  If  aQ  is  chosen  to  be  very  small  (but 
positive)  then  the  5^  order  response  approaches  that  of  the  4^  order 
network  with  k-1  high-pass  elements. 

With  this  understanding  of  the  odd-order  gain  function  characteris¬ 
tics  we  can  again  return  to  the  issue  of  gain-bandwidth  performance.  For 
the  low-pass  case,  where  the  integral  relation  of  eq  (2-11)  applies*  we 
anticipate  that  the  superior  performance  of  the  5th  order  network  in  the 
very  large  region  above  the  passband  (  l  <  x  <  oo  )  will  outweig^  the 
effect  of  the  inferior  performance  below  the  passband  to  result  mo^t  often 
in  the  5^  order  network's  GBW  performance  surpassing  that  of  the 
corresponding  4ta  order  network  in  the  low-pass  constraint  case.  We 
observe*  however*  that  it  is  conceivable  that  for  some  combinations  of 
gain  function  specifications  and  oQ  it  might  be  possible  to  obtain  5th 
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order  GBW  performance  inferior  to  that  for  the  4th  order  network. 

In  the  high-pass  constraint  case,  the  integral  relation  of  eq  (2-12) 

2 

applies,  and  we  notice  that  a  (  1 /u  )  weighting  function  has  now  been 
included.  This  greatly  increases  the  significance  of  the  gain  curve  for 
very  small  values  of  u  (the  region  below  the  passband) .  and  we  observe 
that  this  greatly  increases  the  possibility  of  the  GBW  performance  of  a 
5th  order  network  being  worse  than  that  of  the  4th  order  network  —  for 
the  HP  constraint  case.  In  practice,  we  find  that  this  indeed  occurs, 
with  the  best  GBW  performance  of  a  5t^1  order  network  with  a  HP  constraint 
often  naL  exceeding  that  of  the  corresponding  4ttl  order  network. 

This  discussion  applies  equally  well  to  the  3rd  order  case  as  well  as 
to  all  other  odd-order  networks  synthesized  via  this  technique  of  using  N 
critical  points  in  the  approximation.  In  the  3  order  case,  however,  the 
improvement  in  the  response  above  the  passband  usually  outweighs  the 
poorer  response  below  the  psssband  --  even  in  the  HP  constraint  case  --  so 
that  a  3rd  order  network  optimized  for  gain-bandwidth  will  usually  be  able 
to  exceed  the  GBW  performance  of  a  corresponding  2nd  order  network  regard¬ 
less  of  gain  function  specifications. 

In  summary,  since  the  gain  response  of  an  N  order  odd  network  does 
not  fall  off  as  rapidly  as  that  of  an  N-l  order  even  network  in  the  region 
below  the  passband .  and  since  the  gain -bandwidth  product  is  a  weighted 
integral  over  all  positive  x.  the  order  odd  network  will  in  some  cases 
be  able  to  achieve  gain-bandwidth  performance  superior  to  that  of  a 
corresponding  N-l  order  even  network,  but  in  some  cases  it  will  not.  The 
odd  network  will  always,  however,  be  capable  of  at  least  equalling  the  GBW 


performance  of  the  even  network.  In  addition,  it  is  sometimes  possible  to 
obtain  an  N^  order  odd  network  with  GBW  performance  inferior  to  that  of 
an  N-l  order  even  network  —  this  is  especially  true  in  the  high-pass  con¬ 
straint  case  with  low  gain  slope. 

2.JJL.  Specific  Odd-Order  Gain-Bandwidth  Rfiaultfl 

This  section  presents  the  results  of  specific  studies  into  the 
gain-bandwidth  performance  of  3r<*  and  5^  order  networks  compared  to  that 

of  even-order  networks  with  identical  gain  function  specifications. 

Cases  were  examined  with  0  and  6  dB  per  octave  gain  slopes,  narrow  and 
wide  bandwidths.  and  low-pass  as  well  as  high-pass  reactive  constraints. 
The  plots  illustrate  how  gain-bandwidth  varies  with  the  ripple  specifica¬ 
tion  for  networks  of  different  order. 

The  odd-order  curves  in  all  plots  show  only  the  performance  when  DPC 
is  chosen  to  maximize  gain-bandwidth.  In  all  odd-order  cases  it  is  also 
possible  to  achieve  poorer  GBW  performance  by  choosing  non-optimum  DPC 
values,  if  that  is  desired.  The  GBW  performance  of  the  Nth  order  odd 

networks  can  always  be  reduced  to  that  of  the  corresponding  N-l  order 
even  network  by  appropriate  choice  of  DPC.  and  in  the  aero  gain  slope 
cases  it  is  also  usually  possible  to  obtain  GBW  performance  inferior  to 
th*  N-l  case  as  well. 

Figure  2.13  shows  the  comparative  gain-bandwidth  (GBW)  performance 
of  the  even-  and  odd-order  networks  for  the  low-pass  constraint  case 
(first  reactive  element  is  a  shunt  capacitor).  All  networks  examined  for 


NORMALIZED  GBW  (C, ) 
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Figure  2.13.  Gain-Bandwidth  Results  (GS=0»  BW=1.25.  LP) 

this  case  were  synthesised  for  0  minimum  insertion  loss  (MIL).  0  gain 
slope  (GS).  a  bandwidth  ratio  (BW)  of  1.25*  and  1  high-pass  (HP)  element. 
The  plotted  GBW  data  represents  the  normalised  sise  of  the  first  (shunt 
C)  element  in  the  networks. 

For  this  case  the  odd-order  networks  performed  as  would  be  expected 
from  past  results  with  even-order  networks.  The  GBW  curve  for  the  37^ 

order  lied  between  the  2n<*  and  4**1  order  curves*  and  the  Sth  order  curve 
fell  between  the  A*"*1  and  6t*'  order  curves. 


NORMALIZED  GBW  (l/C, ) 
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Figure  2.14  shows  Che  comparative  GBW  performance  for  Che  high-pass 
constraint  case  (first  reactive  element  is  a  series  capacitor).  All 
these  networks  were  synthesized  for  GS=0,  MIL=0»  BW=1.25.  and  2  HP  ele¬ 
ments.  (The  6^  order  curve  could  not  be  extended  to  higher  ripples 
since  it  was  impossible  to  generate  a  valid  gain  function  that  met  those 
gain  specifications  with  the  specified  topology.) 

The  3rd  order  network  performance  was  better  than  that  of  the  2nd 
order  network,  but  the  4**1  and  St^1  order  networks  exhibited  identical  GBW 


Figure  2.14.  Gain-Bandwidth  Results  (GS=0,  BW=1.25.  HP) 
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performance.  The  s^(p)  functions  derived  from  the  4C**»  5th.  and  6t*1 

order  gain  functions  in  this  case  were  examined  in  detail  to  illustrate 
why  the  5C^  order  GBW  was  approximately  the  same  as  that  for  the  4th 

order  network.  Table  2-1  lists  the  complex  roots  of  the  numerator  and 
denominator  polynomials  of  the  s^(p)  functions  for  each  of  the  three 
networks  in  question  for  the  0.05  dB  ripple  case. 

Comparing  first  the  N=4  and  N=6  root  patterns  (cases  where  there  was 
a  significant  difference  in  GBW  performance),  we  can  see  that  there  are  2 
complex  conjugate  pairs  of  roots  for  both  numerator  and  denominator  in 
the  N=4  case,  while  there  are  3  root  pairs  for  each  in  the  N=6  case. 
Furthermore,  none  of  the  6**1  order  roots  lie  very  close  to  any  of  the  41*1 

Table  2-1.  Roots  of  s^(p)  Function 
RIP  =  0.05  dB  GS  =  0  BW  =  1.25  MIL  =0  HP  Constraint 


Numerator 

Denominator 

N=4 

-0 .002132750 

^0. 82  6  46  94 

-0.1003646 

♦0.7593087 

-0.002497883 

+0.9679665 

-0.1368710 

+1 .0354980 

N=5 

-0.002139014 

+0.8265580 

-0.1008672 

+0.7592448 

-0.002490566 

+0.9680869 

-0.1360134 

+1 .0358420 

-4.557193000 

0.0 

-4.5572740 

0.0 

N=6 

-0.001053324 

+0.8128668 

-0.0459729 

♦0.7  7  97  239 

-0.002188380 

+0.8914190 

-0.0973306 

+0.8964114 

-0.001123381 

+0.9862027 

-0.0503  410 

+1 .0210770 
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order  roots  in  the  complex  plane.  We  also  know  that  the  in-band  gain 
response  of  the  6^  order  network  is  quite  different  from  that  of  the  4ttl 
order  response,  since  the  6**1  order  case  was  approximated  with  7  critical 
points,  and  the  order  case  with  5. 

t’h  th 

In  contrast,  both  the  4  and  5  order  gain  responses  were  approxi¬ 
mated  with  S  critical  points,  and  so  are  almost  identical  within  the 

specified  passband.  In  Table  2-1  we  can  see  that  in  addition  to  the  two 

complex  conjugate  pairs  of  roots  in  the  numerator  and  denominator  of  the 
41*1  order  function,  the  5^  order  function  only  adds  1  real  root  to  each. 

And.  significantly,  the  complex  conjugate  roots  in  the  5^  order  case 

have  been  displaced  only  very  slightly  from  their  positions  in  the  4th 
order  case.  Table  2-II  quantifies  this  almost  imperceptible  change  in 
complex  root  positions  between  the  4th  and  5th  order  caaes. 


Table  2-II.  Percent  Change  in  s^(p)  Root  Positions.  N=4  vs  N=5 

(Table  2-1  Functions) 


Real 

Part 

Imaginary 

Part 

Numerator 

0.294  Z 

0.0108  Z 

0.294 

0.0124 

Denominator 

0.501  Z 

0.0084 

0.631 

0.0332 
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This  represents  an  example  of  the  general  case  discussed  in  the  pre¬ 
vious  section  where  under  the  high-pass  reactive  constraint  the  best  5^ 
order  GBW  performance  is  sometimes  achieved  by  approximating  a  4tn  order 
response • 

Figures  2.15  and  2.16  show  the  GBW  results  for  25Z  bandwidth  net¬ 
works  with  6  dB  per  octave  gain  slope.  We  expect  the  gain-bandwidth  pro¬ 
ducts  of  these  networks  to  be  consistently  greater  than  those  of  the  pre¬ 
vious  discussion  due  to  the  addition  of  gain  sloping,  and  this  is  indeed 
the  case.  Once  again,  the  general  gain- bandwidth  trends  for  the  odd- 
crder  networks  are  as  expected . 

The  gain-bandwidth  results  for  octave  band  designs  with  zero  gain 
slope  are  shown  in  Figure  2.17  (low-pass  constraint)  and  2.18  (high-pass 
constraint).  Once  again*  the  4th  and  5th  order  reaults  in  the  HP  case 

were  identical,  as  they  were  for  the  25Z  bandwidth  case. 

Figures  2.19  and  2.20  show  GBW  results  with  networks  designed  for  6 
dB/ octave  gain  slope  and  octave  bandwidth.  Several  of  the  networks 
exhibited  gain  overshoot  outside  the  specified  passband  for  the  smaller 
ripple  values;  the  GBW  results  for  these  cases  cannot  be  properly  com¬ 
pared  to  the  MIL=0  results  for  the  other  networks,  so  GBW  results  for 
these  overshoot  cases  are  not  included  in  Figures  2.19  and  2.20.  Figure 
2.19  shows  the  results  for  the  low-pass  constraint  case. 

Figure  2.20  shows  the  GS=6.  BW=2.0  results  for  the  high-pass  con¬ 
straint  case.  For  this  case  the  odd-order  networks  performed  as  expected 
in  comparison  to  the  even-order  networks. 
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Figure  2.15.  G»in-Bandwidth  Result*  (GS=6.  BW=1.25,  LP) 
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Figure  2.19.  Gain-Bandwidth  Results  (GS=6.  BW=2,  LP) 
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Figure  2.20.  Gain  Bandwidth  Results  (GS=6,  BW=2.  HP) 
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Another  aspect  of  gain- bandwidth  performance  is  that  the  GBW  product 
of  a  network  should  increase  as  the  gain  reduction  (MIL)  in  the  response 
is  increased.  To  illustrate  this,  the  ripple,  bandwidth,  and  gain  slope 
specifications  are  held  fixed  while  synthesizing  networks  with  a  varying 
MIL  specification  and  noting  the  resulting  GBW  products. 

Figure  2.21  shows  the  GBW  results  as  MIL  is  varied  for  2a<*  through 
6C^  order  networks  with  0.10  dB  ripple,  zero  gain  slope,  and  octave 
bandwidth.  The  low-pass  constraint  case  is  shown  in  Figure  2.21.  and  the 
odd-order  data  plotted  is  that  for  the  optimum  GBW  DPC  values.  The  data 
shows  that  the  performance  of  the  odd-order  networks  follows  the  same 
pattern  as  that  of  the  even-order  networks,  which  is  as  we  would  expect. 

Figure  2.22  shows  the  GBW  results  as  MIL  is  varied  for  the  high-pass 
constraint  case*  with  the  odd-order  data  again  showing  the  optimum  GBW 
DPC  results.  The  odd-order  networks  again  follow  the  same  general  pat¬ 
tern  as  the  even  networks,  and  we  notice  that,  just  as  in  Figure  2.14 
where  ripple  is  being  varied  while  MIL  is  held  fixed,  the  GBW  products  of 
the  4C**  and  5^  order  networks  are  the  same. 

In  summary,  the  gain-bandwidth  performance  of  odd-order  networks 
follows  the  same  pattern  as  that  of  even-order  networks  with  respect  to 
variations  in  the  ripple  and  minimum  insertion  loss  specifications. 

While  an  N**1  order  odd  network  usually  exhibits  gain-bandwidth  perfor¬ 
mance  somewhat  better  than  a  corresponding  N-l  order  even  network,  in 
some  cases  (especially  with  the  high-pass  reactive  constraint  and  low 
gain  slope),  the  odd  network  will,  at  best,  only  equal  the  gain-bandwidth 
performance  of  the  even  network. 
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2.12.  Norton  Transformation 

For  those  instances  where  a  3r<*  order  network  cannot  be  found  to 

exactly  absorb  an  FET's  parasitic  capacitance  while  providing  a  proper 
impedance  transformation.  it  may  be  necessary  to  choose  a  4^  order 

topology  that  allows  a  Norton  Transformation  [8]  to  solve  the  problem. 

When  a  Norton  Transformation  is  to  be  used,  the  strategy  is  to  syn¬ 
thesize  the  matching  network  so  that  it  absorbs  the  FET  parasitic  capaci¬ 
tance  exactly,  but  without  regard  to  the  required  terminating  impedance. 
Then,  a  Norton  Transformation  (Figure  2.23)  is  performed  on  two  of  the 
circuit  elements  to  adjust  the  terminating  impedance  to  50  ohms  (or  some 
other  desired  level). 

While  there  are  occasions  when,  due  to  combinations  of  circuit  ele¬ 
ment  values  and  size  of  the  impedance  adjustment  required,  the  Norton 
Transformation  cannot  provide  the  needed  impedance  transformation,  in 
many  cases  it  provides  a  very  convenient  method  to  adjust  the  output 
impedance  of  a  matching  network. 

Figure  2.23(a)  represents  an  initial  network  with  an  unsatisfactory 
output  impedance  —  here  it  has  been  drawn  with  an  ideal  transformer  to 
bring  the  impedance  up  to  the  desired  R  ohms.  Since  the  network  contains 
adjacent  series  and  shunt  impedances  of  the  same  type  (usually  inductors, 
but  capacitors  would  work  as  well),  we  can  perform  a  Norton  Transforma¬ 
tion  to  eliminate  the  ideal  transformer.  The  procedure  is  to  "split"  2, 


into  two  parts  (Z^  and  Zlb>  such  that 
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where  N  is  the  turns  ratio  of  the  ideal  transformer  to  be  eliminated  from 
the  network.  Then  by  changing  the  topology  to  that  shown  in  Figure 
2.23(b)  where 


and 


Z  =  Z  +  z 

A  lb  2 


Zlb  *  <  Z.b  *  Z2> 


the  transformer  may  be  eliminated  from  the  circuit • 


2>11*  3rd.  Order  re  k tk  Order  Network  Performance 

A  2nd  order  network*  consisting  of  a  single  inductor  connected  at 

the  input  or  output  of  an  FET*  will  almost  never  be  capable  of  simultane¬ 
ously  absorbing  an  FET's  parasitic  capacitance  while  providing  the  proper 
impedance  transformation  to  match  the  device  to  its  load.  In  the  past* 
this  has  meant  that  the  simplest  practical  matching  network  that  could  be 
synthesized  was  of  4**1  order. 

Mow  that  odd-order  networks  can  be  synthesized*  it  is  usually  possi¬ 
ble  to  simultaneously  absorb  FET  parasitics  and  provide  required 
impedance  transformations  with  a  3rd  order  network*  using  just  1  inductor 
and  1  capacitor*  rather  than  using  a  46*1  order  network  requiring  3  induc¬ 
tors  and  a  capacitor.  The  reduction  in  required  chip  area  is  very 
dramatic  for  monolithic  realizations,  and  for  most  applications  the 
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performance  of  the  3  order  networks  will  be  quite  acceptable. 

In  the  examples  that  follow,  the  3  rd  order  networks  were  designed  by 
varying  both  the  ripple  specifications  and  the  denominator  polynomial 
constant  (DPC)  until  a  combination  was  found  that  simultaneously  absorbed 
the  FET  parasitic  and  produced  the  needed  impedance  tranaformation  to 
match  the  FET  to  a  50  ohm  source  and  load.  The  4**1  order  networks  were 
designed  by  varying  the  ripple  and  bandwidth  specifications  until  a  com¬ 
bination  was  found  that  absorbed  the  FET  parasitica.  and  then  applying  a 
Morton  Transformation  to  the  inductors  to  adjust  the  terminal  impedance 
to  50  ohms.  (In  one  case  the  Morton  Transformation  was  not  necessary.) 


2.11.1.  Low-Mo i a Amplifier  Example 

This  section  will  compare  the  performance  of  3rd  and  4th  order 
matching  networks  used  in  low-noise  amplifiers  (LNA's)  by  presenting 
actual  input  and  output  matching  network  designs  for  a  10.4  -  12.4  GHz 
LNA  using  the  Cornell  513-4  0.5  *im  low-noise  FET. 

Figure  2.24(a)  shows  the  input  model  of  the  device,  optimized  for 
best  noise  match  over  10.4  -  12.4  GHz.  The  matching  networks  shown  in 
Figure  1  both  exactly  match  this  model  to  50  ohms.  Both  networks  were 
synthesized  for  1  dB  per  octave  gain  slope  (GS)  and  a  minimum  insertion 
loss  (MIL)  of  0  dB.  While  the  4C^  order  (N=4)  design  provides  a  good 

match  over  a  wider  bandwidth  (BW)  than  does  the  3rd  order  network,  this 
extra  bandwidth  ia  not  needed  in  this  application.  The  4C^  order  design 
also,  however,  requires  2  more  inductors  than  does  the  3rd  order  network. 


(a) 
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Figure  2.24.  Low-Noise  Amplifier  Input  Networks 


and  one  of  theae  ia  quite  large  for  12  GHz  operation  (1.8  nH) . 

Figure  2.25(a)  compare*  the  inaertion  loaa  of  the  two  LNA  input  net¬ 
work*  (analyzed  with  the  apeeified  model),  and  Figure  2.25(b)  eomparea 
the  input  noiae  figure  obtained  with  the  two  networke.  Theae  reaulta 
ahow  that  the  performance  of  the  3r^  order  network  ia  quite  comparable  to 
that  of  the  4th  order  network  over  thia  17Z  band,  and  that  at  ita  worat 
point  the  3rd  order  deaign  only  differs  from  the  A11*1  order  design  by 
about  0.1  dB  in  both  insertion  loss  and  noise  figure. 

Figure  2.26(a)  shows  the  output  model  of  the  513-4  device  (with  the 
network  of  Figure  2.24(c)  connected  at  its  input).  Again,  both  output 
networks  in  Figure  3  were  synthesized  to  match  this  model  to  50  ohms. 

Once  again,  the  4^  order  network  provides  a  match  over  a  wider  bandwidth 
than  needed,  and  requires  two  more  (fairly  large)  inductors  than  does  the 
3rd  order  design. 

Figure  2.27  compares  the  insertion  loss  performance  of  the  two  out¬ 
put  networks,  and  the  3r<*  order  network  again  performs  comparably  with 
the  4th  order  network,  but  with  somewhat  greater  ripple  apparent  in  the 
response. 

2.11.2.  Power  Amplifier  Rrampie 

In  general,  the  matching  problems  presented  by  a  power  FET  are  quite 
different  from  those  of  a  low-noise  FET.  due  to  the  very  different  resis¬ 
tance  and  reactance  values  required  to  model  these  physically  very  dif¬ 
ferent  devices.  For  completeness,  this  section  will  compare  the  perfor- 
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Figure  2.25.  Low-Noise  Amplifier  Input  Network  Performance 
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Figure  2.26.  Low-Noise  Amplifier  Output  Networks 
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Figure  2.27.  Low-Noise  Amplifier  Output  Network  Performance 


0.605  pF 


POWER  FET 
INPUT  MOOEL 
8.5 -9. 5  GHz 

(CORNELL  1505  0EVICE)O 


7.45  ft 


Figure  2.28.  Power  FET  Input  Model 


•x 


42 


nance  of  3r<*  and  4^  order  input  and  output  networks  when  designed  for 
use  with  a  high-power  FET  amplifier  (using  an  output  model  for  the  device 
obtained  under  large  signal  conditions)* 

Figure  2.28  shows  the  input  model  of  the  Cornell  1505  0.8  pm  power 
FET*  optimized  over  8.5  -  9.5  GHz  for  best  match  to  the  measured  RF  data. 
Both  input  networks  of  Figure  2.29  match  this  model  to  a  50  ohm  source 
with  6  dB  per  octave  gain  slope.  All  elements  of  the  3rd  order  network 
are  realizable*  but  the  4th  order  network  at  (b)  requires  an  8.9  nH 
inductor,  which  is  much  too  large  to  fabricate  at  9.5  GHz.  The  only 
other  4th  order  topology  that  can  potentially  allow  easy  bias  introduc¬ 
tion  (without  adding  a  large  RF  choke)  is  shown  at  (c).  The  response  of 
this  network  is  essentially  the  same  as  that  for  the  network  at  (b)»  but 
the  shunt  inductor  in  (c)  (0.045  nH)  is  too  small  to  allow  direct  intro¬ 
duction  of  gate  bias.  Therefore,  and  for  different  reasons*  neither  of 
the  4th  order  input  topologies  provide  a  practical  monolithic  solution  to 
our  matching  problem. 

Figure  2.30  compares  the  gain  responses  of  the  3r<*  and  4**1  order 
input  networks.  While*  in  this  case*  the  4th  order  response  is  closer  to 
an  "ideal"  sloped  response  than  is  that  of  the  3r<*  order  network*  it  is 
not  realizable  monolithically*  and  the  3r<*  order  response  would  be  a 
satisfactory  starting  point  for  optimization  in  most  applications. 

The  large  signal  output  model  for  the  device  (at  an  output  power 
level  of  +21  dBm)  is  shown  in  Figure  2.31*  and  the  output  matching  net¬ 
works  are  shown  in  Figure  2.32.  Both  networks  were  designed  for  a  flat 
sutch  across  the  band*  and  the  ripple  specifications  are  very  similar. 
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Figure  2.29.  Power  Amplifier  Input  Networks 


FREQ  (GHz) 


0.475  pF 


LARGE  SIGNAL 
OUTPUT  MOOEL 
8.5 -9.5  GHz 
(CORNELL  1505  DEVICE) 


Figure  2.31.  Power  FET  Output  Model 
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Figure  2.32.  Power  Amplifier  Output  Network* 


Mo  Norton  Transformation  was  necessary  for  this  4^  order  network*  so 
only  1  additional  inductor  is  required. 

Figure  2.33  shows  the  gain  performance  of  both  networks;  there  is 
no  significant  difference  in  the  two  responses  (less  than  0.03  dB*  worst 
case*  across  the  band). 

2.11.2.  Summary 

Very  often*  a  narrowband  FET  matching  problem  can  be  solved  with  a 
3r<*  order  network  instead  of  with  one  of  4^  order.  Although  this  is  not 


FREQ  (GHz) 


Figure  2.33.  Power  Amplifier  Output  Network  Performance 


always  possible*  the  3rc*  order  networks*  when  usable*  provide  matching 
networks  requiring  substantially  less  chip  area  than  corresponding  4th 
order  networks  with  only  slight  degradation  in  circuit  performance. 

In  practical  monolithic  GaAs  amplifiers*  the  smaller  sire  of  the  3r<* 
order  network*  when  compared  to  a  corresponding  4th  order  network,  will 
almost  always  greatly  overshadow  any  slight  performance  penalty  incurred. 


Chapter  3 


Distributed  Element  Broadband  Ictvork  Design 
(Arbitrary  Order) 


In  some  monolithic  applications  distributed  element  matching  networks 
are  more  convenient  to  use  than  lumped  element  networks.  For  these 
instances  it  would  be  advantageous  to  be  able  to  use  odd-order  networks, 
just  as  it  is  with  lumped  elements.  Our  goal  is  still  to  design  broadband 
matching  networks  with  the  minimum  number  of  elements,  and  once  the  theory 
has  been  extended  to  allow  odd-order  networks  this  will  often  mean  3rd 
order  networks . 

This  chapter  will  describe  the  necessary  modifications  to  the  lumped 
element  synthesis  procedure  to  permit  synthesis  of  odd-  or  even-order  dis¬ 
tributed  element  networks.  These  results  represent  an  extension  of  the 
existing  techniques  which  are  limited  to  even-order. 


l.L.  Equiripple  Cain  Enaction  Synthes ia  -  Even  qx.  DM  Order 

With  a  few  exceptions,  the  overall  network  synthesis  algorithm 
described  in  g  2.1  for  lumped  elements  remains  valid  for  distributed  ele¬ 
ments.  Rather  than  conducting  the  gain  function  synthesis  in  the  fre- 
quency  domain  (u  )•  however,  for  distributed  element  networks  it  is  con- 

ducted  in  the  Richards1  Transformation  domain  (  Q  )  as  discussed  below. 
The  other  change  is  that  the  element  extraction  step  is  conducted  dif- 
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ferently  in  order  Co  realize  the  network  as  shorted  and  open-circuited 
stubs  and  transmission  line  sections  instead  of  as  lumped  inductors  and 
capacitors . 

The  distributed  element  networks  to  be  synthesized  are  assumed  to 
consist  of  an  arbitrary  cascade  of  commensurate  length  unit  elements  (cas¬ 
cade  transmission  line  sections)  as  well  as  series  or  shunt  open-circuited 
or  short-circuited  stubs.  The  total  number  of  circuit  elements  may  be 
even  or  odd. 

Gain  function  synthesis  is  carried  out  in  the  Richard's  Transforma¬ 
tion  domain  C  93 


Q  =  tan  wT 


(3-1) 


where  u  is  frequency •  and  T  is  the  delay  length  of  the  commensurate  lines. 
The  transfer  function  of  the  distributed  element  network  can  be  con¬ 
veniently  expressed  in  the  following  form  [10] 


G(x) 


K*kCl  +-x),.q 

Vx) 


(3-2) 


o 

where:  x  =  ft  .  the  square  of  the  Richards'  Transformation  variable 

K  =  gain  constant  chosen  to  maintain  realizability 
(0  5  G(x)  SI  V  x  >  0  ) 

k  =  number  of  high-pass  elements  in  the  network 

(known  once  the  network  topology  has  been  specified) 

q  -  number  of  unit  elements  in  the  network 

(known  once  the  network  topology  has  been  specified) 

P„(x)  =  an  Nth  order  polynomial  in  x »  whose  coefficients 
are  to  be  determined 
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This  expression  is  slightly  different  from  the  G(x)  form  used  with 

lumped  elements  (eq  (2-1))  •  but  a  similar  procedure  will  be  used  to  form 

an  equiripple  approximation  to  a  specified  gain  function  by  determining  K 

and  the  coefficients  of  P„(x). 

N 

1.1.1.  Approximation  Algorithm 

Since  the  general  gain  function  approximation  strategy  for  distri¬ 
buted  elements  is  the  same  as  that  for  lumped  elements*  the  detailed 
explanation  of  the  algorithm  presented  in  g  2.4.2  will  not  be  repeated 
here.  This  section  will  primarily  summarize  the  modified  equations  used 
in  the  computations. 

As  in  the  lumped  element  case*  the  approximation  is  specified  over  N 
critical  points  when  N  is  odd*  and  N+l  critical  points  when  N  is  even. 
After  choosing  the  initial  estimate  of  the  critical  points  (x^) .  setting 

K=l.  and  computing  the  required  gain  levels  at  each  critical  point*  we 
have 
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For  odd  N,  we  again  choose  Bq  as  a  parameter,  and  eq  (3-4)  then 
represents  a  system  of  N  equations  in  N  unknowns*  since  ()  i-8  speci¬ 
fied  at  N  critical  points.  For  even  N,  an  is  also  treated  as  a  variable; 
the  N+l  critical  points  specified  allow  determination  of  all  N+l  coeffi¬ 
cients  of  PN(x),  and  eq  (3-4)  represents  N+l  equations  in  N+l  unknowns. 

Equation  (3-4)  may  be  conveniently  represented  in  matrix  notation  by 


X  a  =  b 


where,  for  odd  ft. 


X 


N  N-l  2 

^  XN  *  XN  *N  ] 


(3-5) 


(3-6) 


aN*  aN-l . 


V  ai 


1T 

J 


and 


(3-7) 


_  j  .^(1^)1 

'  J  Wl>  0  •  <W*2>  ° 


•  Gr.,<V  0  J 


(3-8) 
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For  even  ve  have 


f  N+l  N 

X1  X2 

N+l  N 


X  = 


N+l 

XN 


N 


N+l  N 


X1  1 


x2  1 


X  1 

n 


L  ’Sl+l  XN+1  *  ^N+l 


(3-9) 


fl  =  |  a^ •  a^»  Aq  J 


(3-10) 


:i^(l+x1)q 

.(*,) 


req 


x^(l+x7)q 

Greq<X2> 


4(uVq  -iS«i(1<Vi1'1  j  (3-u) 

Greq^XN^  *  Greq(xN+l*  j 


We  again  solve  eq  (3-5)  by  LU  factorization  for  best  numerical  sta¬ 
bility. 


This  completes  the  first  part  of  the  approximation  algorithm,  and  the 
coefficient  vector,  a.  is  now  available.  The  remainder  of  the  procedure 
deals  with  determining  the  new  critical  points  of  the  response  via  the 
Newton-Raphson  Method,  and  this  is  accomplished  in  exactly  the  same  manner 
as  it  is  for  lumped  elements  as  described  in  g  2. A. 2  (only  the  functions 
and  derivatives  evaluated  during  the  procedure  are  different). 


1.1.2,.  Implementation 

The  above  algorithm  has  been  implemented  in  a  computer  subroutine  for 
solving  the  distributed  element  approximation  problem  with  functions  of 


arbitrary  order  and  combined  with  other  existing  subroutines  to  form  a 
complete  distributed  element  network  synthesis  computer  program  (desig¬ 
nated  AD 1ST) .  A  commented  FORTRAN  listing  of  the  approximation  subroutine 
(designated  ODDIST)  is  presented  in  Appendix  A. 


1.2.  Effects  af.  Parying  the  Denominator.  Polynomial  Constant  (QE£) 


The  general  remarks  made  in  g  2 .5  in  reference  to  the  realizability 
constraints  on  the  choice  of  DPC  in  the  lumped  element  case  are  also  valid 
for  the  distributed  element  case.  DPC  must  always  be  chosen  to  be  posi¬ 
tive.  and  there  will  be  some  maximum  DPC  value*  dependent  upon  the  partic¬ 
ular  gain  function  specifications  in  each  case,  which  results  in  a  realiz¬ 
able  gain  function.  (The  results  presented  in  Chapter  4  show  how  this 
maximum  realizable  DPC  value  may  be  determined  in  advance  for  3r<*  order 
networks .) 

The  following  two  sections  present  some  examples  of  the  effects  of 
DPC  variation  in  actual  distributed  element  networks. 

2*2.1.  Gain-Bandwidth  BgjuULioi  with  Q££ 

Odd-order  distributed  element  networks  exhibit  gain-bandwidth  (GBV) 
variations  as  the  value  of  DPC  is  changed,  just  as  described  for  the 
lumped  element  networks.  The  exact  nature  of  this  variation  is  very 
dependent  upon  the  combination  of  gain  function  specifications  for  each 
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particular  case*  so  no  general  statement  can  be  made  as  to  the  overall 
shape  of  the  GBW  vs  DPC  curves. 

Figure  3.1(a)  shows  an  example  of  GBW  behavior  for  the  low-pass  con¬ 
straint  case*  and  Figure  3.1(b)  shows  an  example  of  the  results  for  the 
high-pass  constraint  case.  Cases  also  exist  (not  shown)  where  GBW  reaches 
a  maximum  at  some  DPC  value  internal  to  the  realizable  interval*  just  as 
with  lumped  elements;  the  two  cases  shown  in  Figure  3.1  reach  maximum  GBW 
at  one  or  the  other  end  of  the  realizable  interval.  Both  examples  in  Fig¬ 
ure  3.1  are  3rd  order  networks  synthesized  for  octave  bandwidth*  0.1  dB 
ripple*  zero  gain  slope*  and  zero  minimum  insertion  loss. 


1.2.2.  Terminating  Impedance  Behavior  with  DPC 

Just  as  with  lumped  elements*  changing  DPC  in  an  odd-order  distri¬ 
buted  element  network  also  alters  the  required  terminating  resistance  for 
the  network.  Figure  3.2  shows  two  examples  of  this  behavior*  for  the  same 
two  circuits  (and  same  gain  specifications)  used  to  generate  the  gain 
bandwidth  curves  in  Figure  3.1.  Once  again*  the  only  general  statement 
which  may  be  made  regarding  required  terminating  resistance  as  DPC  is 
changed  is  that  it  varies*  the  exact  shape  of  the  curve  is  dependent  upon 
the  particular  combination  of  gain  specifications  in  each  case* 


1.1.  Chpoaiog  Q££ 

The  denominator  polynomial  constant  is  used  as  an  additional  syn¬ 
thesis  parameter  when  designing  odd-order  distributed  element  matching 
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1-A.  Eraciical  Monolithic  Input  networks 

Usable  monolithic  topologies  for  distributed  element  networks  are 
rather  limited  *  since  we  can  only  allow  shunt  open-  or  short-circuited 
stubs  and  cascade  transmission  line  sections.  Series  open-  or  short- 
circuited  stubs  cannot  be  realized  in  microstrip*  and  are  therefore  not 
included  in  any  circuits  intended  for  monolithic  applications.  (Series 
open-circuited  stubs  are  included  in  some  of  the  topologies  discussed*  but 
only  to  model  the  FET  input  reactance.  These  elements  are  not  physically 
constructed .) 

We  require  the  circuits  to  include  a  shunt  shorted  stub;  forming  the 
"short"  with  a  bypass  capacitor  will  then  produce  a  convenient  point  for 
bias  introduction.  Distributed  synthesis  does  not  allow  including  series 
capacitors,  so  a  DC  blocking  capacitor  will  need  to  be  added  to  all  the 
example  networks. 

Figure  3.3  shows  several  example  distributed  input  network  topologies 
that  can  be  realized  monolithically.  There  are  other  possible  topologies* 
as  well*  but  these  represent  the  simplest  ones  that  fulfill  our  basic 
requirements.  The  series  distributed  FET  input  model  is  shown  at  the 
right  of  each  circuit  in  the  figure*  and  the  assumed  resistive  source 
impedance  for  the  amplifier  is  shown  at  the  left. 

The  3r<*  order  network  shown  in  Figure  3.3(a)  is  the  simplest  topology 
which  is  likely  to  be  able  to  simultaneously  absorb  the  FET  parasitic 
while  producing  the  correct  impedance  transformation.  Figures  3.3(h) 
and  (c)  are  5**1  order  networks*  and  the  only  difference  between  these  two 
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for  a  fixed  gain  function  specification  will  be  the  values  of  the  required 
terminating  impedance.  Figures  3.3(d)  and  (e)  are  4^  order  networks. 

Several  of  the  networks  in  Figure  3.3  also  have  the  potential  of 
allowing  a  Kuroda  Transformation  (discussed  in  g  3.8)  to  further  adjust 
the  terminating  resistance. 


1.1.  Practical  Monolithic  Output,  networks 

Figure  3.4  shows  several  example  distributed  output  network  topolo¬ 
gies  that  can  be  realized  monolithically •  The  same  basic  circuit  con¬ 
straints  discussed  in  g  3.4  also  apply  to  these  networks.  The  shunt  dis¬ 
tributed  FET  output  model  is  shown  at  the  left  of  each  circuit  in  the  fig¬ 
ure.  and  the  assumed  resistive  load  impedance  for  the  amplifiers  is  shown 
at  the  right.  While  there  are  other  possible  outpout  topologies*  the  ones 
shown  are  the  simplest  that  fulfill  our  basic  requirements. 

The  circuit  shown  at  Figure  3.4(a)  is  a  3r<*  order  network,  those  at 
(b)  and  (c)  are  4^  order  networks,  and  that  at  (d)  is  5^  order.  All 
these  circuits  have  the  potential  of  allowing  a  Kuroda  Transformation  to 
adjust  the  terminating  impedance. 


inilHil'Hi  -i  1  ilia 
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!•£..  Eractical  Monolithic  Interstage  Heivorfca. 

Figure  3.5  shows  several  examples  of  distributed  interstage  network 
topologies  that  can  be  realized  monolithieally .  The  distributed  FET  out¬ 
put  model  is  shown  at  the  left  of  each  circuit  in  the  figure,  and  the  dis¬ 
tributed  FET  input  model  is  shown  at  the  right. 

The  circuits  shown  at  Figure  3.5(a-c)  are  4  and  5L  order  networks 

that  may  be  directly  synthesized.  The  first  strategy  to  be  used  in 
designing  circuits  of  this  type  is  to  attempt  to  find  a  combination  of 
gain  function  parameters  that  simultaneously  provide  an  impedance 
transformation  from  to  Rj  while  exactly  absorbing  and  Z^.  If  this 
cannot  be  done,  it  may  be  possible  to  solve  the  problem  by  absorbing  Z^^ 

exactly  and  then  using  a  Kuroda  Transformation  to  adjust  the  impedance 
transformation  ratio.  (Z^  will  be  properly  absorbed  in  this  process  if 

the  ratio  of  R^  to  Z^  prior  to  performing  the  transformation  is  the  same 
as  that  required  by  the  FET  input  model.  Since  R2^Z02  *8  i-nvariant  under 
impedance  scaling,  once  the  output  impedance  of  the  network  is  adjusted  to 
Rj  by  the  transformation.  Z02  will  also  be  scaled  to  the  proper  level  for 
exact  absorption.) 

If  neither  of  these  strategies  is  successful,  it  is  often  possible  to 
solve  the  interstage  problem  by  cascading  a  3rd  or  4th  order  network  at 

the  output  of  the  first  FET  with  a  2nd  order  network  at  the  input  to  the 
second  FET.  The  resulting  topologies  are  shown  at  Figure  3.5(d)  and  (e). 
The  strategy  is  the  same  as  with  lumped  elements;  a  2nd  order  network  is 
synthesized  st  the  input  to  the  second  FET  by  adding  shunt  shorted  stub  Z* 
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such  that  the  FET  parasitic  is  absorbed  exact ly.  Then  a  3r<*  order  output 

network  is  synthesized  for  the  first  FET  to  match  to  the  impedance  level 
required  at  the  input  to  the  2n<1  order  network.  Due  to  the  loading  and 
ripple  effects  between  the  two  networks,  it  may  be  necessary  to  apply  some 
computer  optimization  to  the  cascaded  circuit  to  obtain  a  satisfactory 
overall  gain  response,  especially  for  broadband  networks. 

As  with  lumped  element  interstage  networks,  it  is  again  important 
that  all  possible  distributions  of  the  zeroes  of  the  8j^(p)  function  be 

examined  during  the  search  for  a  suitable  interstage  network. 

1.1.  Gain-Bandwidth  Results 

Odd-order  distributed  networks  exhibit  gain-bandwidth  (GBW)  charac¬ 
teristics  similar  to  those  of  lumped  element  networks.  The  detailed  dis¬ 
cussion  of  8  2.11  will  therefore  not  be  repeated,  but  two  example  plots 
will  be  presented. 

Figure  3.6  shows  the  normalized  GBW  product  of  3r<*»  4**1.  5^.  and  6t^1 
order  distributed  networks  plotted  against  the  ripple  specification  for 
the  low-pass  constraint  case.  The  GBW  performance  of  the  3  and  5^ 
order  networks  is  much  as  we  would  expect  compared  to  the  4^  and  6th 
order  networks,  and  is  very  similar  to  that  shown  in  Figure  2.17  for 
lumped  elements. 

Figure  3.7  shows  the  GBW  results  for  a  set  of  distributed  networks 
for  the  high-pass  constraint  case.  Again,  the  results  are  very  similar  to 
the  lumped  element  high-pass  constraint  results  shown  in  Figure  2.18.  The 
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Figure  3.6  Gain-Bandwidth  Results  (LP»  Distributed) 


4th  and  5th  order  GBW  curves  overlie  each  other,  again  because  the  root 
patterns  of  the  s^Cp)  functions  are  so  nearly  identical.  (The  6th  order 

curve  could  not  be  extended  to  higher  ripple  values  since  it  was  impossi¬ 
ble  to  generate  a  valid  gain  function  that  met  those  gain  specifications 
with  the  specified  topology.) 

The  odd-order  results  here  for  both  the  low-pass  and  high-pass  con¬ 
straint  cases  reflect  the  maximum  GBW  possible  for  each  value  of  ripple, 
obtained  by  optimum  choice  of  DPC. 


NORMALIZED  GBW  (Z0I) 
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Figure  3.7.  Gain-Bandwidth  Results  (HP.  Distributed) 


GBW  results  are  not  presented  here  for  odd-order  gain  sloped  distri¬ 
buted  networks  nor  for  the  case  of  fixed  ripple  and  varying  minimum  inser¬ 
tion  loss •  The  distributed  network  curves  are  generally  similar  to  the 
lumped  element  plots  presented  in  Figures  2.19  through  2.22. 

1«&.  Kuroda  Transformation 

For  those  instances  where  a  synthesized  network  cannot  be  found  to 


exactly  absorb  an  FET's  parasitic  reactance  while  providing  a  proper 
impedance  transformation,  it  may  be  possible  to  employ  a  high-pass  Kuroda 
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Transformation  [11]  to  adjust  the  impedance  transformation  ratio  to  the 
required  value • 

Figure  3.8  illustrates  the  high-pass  Kuroda  Transformation.  Figure 
3.8(a)  represents  an  initial  distributed  network  with  an  unsatisfactory 
output  impedance  —  here  it  is  dram  with  an  ideal  transformer  to  bring 
the  output  impedance  up  to  the  desired  R  ohms.  Since  the  network  contains 
an  adjacent  shunt  short-circuited  stub  and  cascade  line  section,  we  may 
employ  a  Kuroda  Transformation  to  eliminate  the  ideal  transformer  from  the 
network. 

The  procedure  is  to  "split"  shunt  shorted  stub  into  two  parts  (zja 
and  Zj^)  such  that 

*2. 

N  =  1  +  7^- 
lb 

where  N  is  the  turns  ratio  of  the  ideal  transformer  to  be  eliminated  from 
the  network.  Then  by  changing  the  topology  to  that  shown  in  Figure  3.8(b) 
where 


and 


Zlb  *  Z2 
Zlb  *  Z2 


_ lib _ 

Zlb  +  Z2 


the  ideal  transformer  may  be  eliminated  from  the  circuit. 
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Chapter  4 


Amalytical  Detcraiaatioa  of  Maiian  Denominator  Polynomial 
Constant  for  Kealinability 


Chapter  2  discussed  the  limited  range  of  the  denominator  polynomial 
constant  (DPC)  that  can  result  in  realizable  network  gain  functions  for 
odd-order  networks.  This  realizable  range  is  bounded  by  zero  and  some 
positive  value  dependent  upon  the  particular  gain  function  specifications 
for  each  case.  It  would  be  desirable  to  be  able  to  determine,  in  advance, 
this  positive  limit  for  DPC  so  the  circuit  designer  would  have  this  infor¬ 
mation  available  at  the  time  of  synthesis. 

This  chapter  presents  techniques  for  analytically  determining  the 
maximum  value  for  DPC  that  produces  a  realizable  gain  function  for  both 
lumped  and  distributed  element  3rd  order  networks  with  arbitrary  topology, 
bandwidth,  ripple,  and  gain  slope.  The  procedure  is  completely  general, 
and  applies  to  all  3rd  order  networks. 


Rli.TBi.nt-n 

Figure  4.1  shows  a  generic  3rd  order  bandpass  gain  function,  of  arbi¬ 
trary  bandwidth,  ripple,  and  gain  slope,  that  illustrates  the  mathematical 
problem  being  solved.  The  three  critical  frequencies  of  the  response  are 
designated  Xj,  x2>  and  x^;  for  the  lumped  element  case  we  can  simplify 
matters  somewhat  by  normalizing  the  problem  so  that  x^l •  Let  the 
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specified  gain  levels  at  these  critical  points  be  ,  G^.  and  G^.  respec¬ 
tively. 

As  in  Chapter  2«  the  3rt*  order  lumped-element  bandpass  gain  function 
to  be  used  is 


G(x) 


Ax  +  Bx  +  Cx  +  D 


(4-1) 


where  the  numerator  gain  constant  (K)  has  been  set  to  1  and  k  is  the 
number  of  high-pass  elements  in  the  network,  (k  =  1  or  2  for  the  bandpass 
networks  of  interest.)  We  wish  to  set  up  a  system  of  equations  using  eq 
(4-1)  and  the  desired  gain  response  specifications  that  allow  us  to 


to 

determine  the  maximum  value  of  D*  above*  that  reaults  in  a  positive  coef¬ 
ficient  for  the  high  order  denominator  term  (A). 

Using  the  specified  gain  levels  at  the  three  critical  points*  we  can 
immediately  write  the  following  three  relations: 


Ax*  +  BXj  +  CXj  +  D 


h 

G, 


(4-2) 


Ax*  +  Bx*  +  Cx2  +  D 


(4-3) 


A  +  B  +  C  +  D  = 


(4-4) 


Uhere  we  have  used  the  fact  that  x^=l  for  this  case. 

A  fourth  equation  is  provided  by  the  constraint  that  the  alope  of  the 
gain  function  at  the  center  critical  point  be  equal  to  the  gain  slope 
specified  for  the  response.  If  we  let  S  be  this  specified  slope*  then  we 
require  that  Gf(x)[c=x  =  s*  It  can  be  shown  that 


rf,„ x  _  (A^W+Cs+P)!^"1  -  xk(3Ax2+ZBx+C) 

®  'X/  “19  5 

(AxJ  ♦  Bx*  +  Cx  +  D)z 


(4-5) 


Evaluating  at  x=*2 •  this  relation  can  be  simplified  considerably  by  using 
eq  (4-3) 


G  g 

lx,,  *  -b  >  -  -b 3A*!  *  2*‘2  * c  > 

2  x2  x2 


(4-6) 


Setting  eq  (4-6)  equal  to  our  specified  slope*  S*  and  re-arranging*  pro¬ 
duces  an  equation  which  may  be  added  to  the  three  above  to  yield  a  solv¬ 
able  system. 


3x‘A  +  2x2B  ♦  C  = 


(4-7) 


By  using  x2  «s  s  paramatar  to  be  chosen  (instead  of  as  a  variable  to 

be  determined) •  we  have  a  linear  system  of  four  equations  in  four  unknowns 
(the  coefficients  A.  B,  C.  D)  which  is  solvable.  This  system  can  be  con¬ 
veniently  represented  in  matrix  notation  by 

r  i 

1*1 


1 

1 

3 

X1 

X1 

3 

X2 

X2 

3x22 

2x 

xi  1 


x2  1 


A  ' 

B 

C 

D  . 


(4-8) 


corresponding  to 


(4-9) 


where 


(4-10) 


Once  x2  is  chosen,  the  matrix  X  is  determined,  and  we  can  calculate 

the  values  for  each  element  of  b  as  soon  as  the  gain  function  is  speci¬ 
fied*  so  we  can  solve  for  the  coefficient  vector,  a.  This  is  not  quite 
the  way  we  will  proceed,  however,  since  we  desire  a  different  type  of 
"solution" . 
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By  using  determinants  we  can  derive  an  analytical  expression  for  the 
value  of  A.  We  can  then  find  the  value  of  x 2  that  makes  A  =  0.  We  know 

that  this  X2  will  be  a  unique  solution  because  the  equiripple  solution  to 

the  approximation  problem  is  itself  unique  as  shown  by  Petersen  Cl]. 
Deriving  the  analytical  expression  for  D*  and  substituting  the  x2  value 

found  to  produce  a  zero  value  for  A  will*  in  turn*  show  the  corresponding 
D  value  that  forces  A  to  zero.  This  value  for  D  (DPC)  then  represents  the 
maximum  value  that  produces  a  realizable  gain  function;  any  larger  DPC 
value  will  produce  a  negative  A  value*  which  makes  the  gain  function 
unrealizable. 


Taking  only  the  numerator  of  the  full  expression  for  A.  and  setting 
it  equal  to  zero*  we  have: 

0  -  J  .Uj  - 1)  ♦  £  -  f  |  4 


L 


'1 


U 


L  x*q  -  V  £  1 

lG„  G„  G, 

+  o(l  -  x\) 

l  3  2  1 

k,  2  k  2  1 

xAx.  -  1)  x.  x. 


!•<*.  C2  *0,  Gj  I 


(4-11) 


where  a  is  the  function  of  S*  x^.  and  defined  in  eq  (4-10). 


Mote  that  since  we  specify  the  reference  gain  function  as  an  exponen¬ 
tial  function  of  x*  the  slope  vslue*  S*  is  also  really  an  exponential 
function  of  x 2.  In  the  general*  arbitrary  gain  slope*  case  G2  is  also  an 

exponential  function  of  x2«  These  complications  make  eq  (4-11)  extremely 
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difficult*  if  not  impossible*  to  solve  using  strictly  analytical  tech¬ 
niques*  but  not  particularly  difficult  to  solve  via  numerical  equation 
solving  techniques  on  a  computer.  So  we  use  a  computer  algorithm  to  solve 
eq  (4-11)  for  that  value  of  x2  (x^  <  <  1  )  which  produces  A  =  0. 

The  full  expression  for  DFC  is  also  quite  complex.  If  we  let 


Num(x2) 

DenCxj) 


[  k  1 

I  x,  x,  I 


Num(x2) 


I  G3  Gj  I  2 


[  k  2  1  1 


3x*(x?  -  x,)  x?  x^  | 


♦  |  .(.j  -  «J)  ♦  | 

L  t  1  3  J  j 


I 

♦  |  a(x3  -  Xj^)  + 

L 

r 

+  |  o(x3  -  x3)  + 


+  G3  "  G1  j  2 


2x^(x,  -  x3)  | 


k,  3  2. 

x_(x,  -  x, ) 


where  a  is  again  as  defined  in  eq  (4-10)*  and 


Den(x2>  =  (xj-l)x2  +  2(1-Xj)x3  ♦  (x3-3xj+3xj-l)x2 


+  2xj(1-Xj)x2  +  x^(xj-l) 


(4-12) 


(4-13) 


(4-14) 
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After  substituting  the  known  values  into  equations  (4-13)  and  (4-14)  the 
resulting  maximum  realizable  value  for  D  can  be  computed  via  eq  (4-12). 


This  procedure  has  been  coded  into  a  set  of  computer  subroutines 
(listed  in  Appendix  A),  and  incorporated  into  the  ARBSYN  computer  program 
where  the  maximum  allowable  DPC  value  for  3rd  order  networks  is  computed 
and  presented  to  the  operator  at  the  appropriate  time. 


4.2.  Distributed  Elements 

The  approach  to  finding  the  maximum  realizable  DPC  value  for  distri¬ 
buted  element  networks  is  the  same  as  that  for  lumped  elements,  except 
that  the  equations  are  slightly  different  and  more  complicated.  Figure 
4.1  is  also  valid  for  the  generic  3rd  order  gain  response  in  the  distri- 
buted  case,  with  the  important  exception  that  we  can  no  longer  simply  nor¬ 
malize  so  that  x^=l.  (With  distributed  elements,  due  to  the  effects  of 

chosen  line  length  and  the  Richards  Transformation,  x^  will  often  be  some 
value  other  than  1.) 

As  in  Chapter  3.  the  3rd  order  distributed  element  bandpass  gain 
function  is 

G(x)  =  -  (4-15) 

AxJ  +  Bx^  +  Cx  +  D 

where  k  is  the  number  of  high-pass  elements  and  q  is  the  number  of  cascade 
lines  in  the  network. 
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The  three  equations  derived  from  the  gain  level  specifications  are 


3  2 

Axj,  +  Bxj  +  Cxj  ♦  D 


x*  (1  +  x,)q 
G, 


3  2 

Ax2  +  Bx2  +  Cx2  +  D 


x*  (1  +  x?)q 


3  2 

Ax^  +  bxj  +  Cxj  +  D 


x*  (1  +  x^)q 


The  derivative  expression  evaluated  at  x2  is 


G'(x) 


kG,  qG, 


x=x2_  x2  (l+x2) 


1 

— - * - -<  3Ax,  +  2Bx,  +  C  ) 

x*(l+x„)q  2  2 


After  setting  eq  (4-19)  equal  to  the  specified  gain  slope*  S*  and 
arranging  terms  we  have 


,  kx^”*( 1  +  x  )q  +  x^q(l  +  x  )q_1 

3Ax2  +  2Bx2  +  C  =  — ‘ - * - - - A 4 - 


-  -S-  f  xk(  l  +  x  )q 

L  *2u  ♦  v  J 


We  again  treat  x2  as  a  parameter*  and  the  resulting  system  of 


(4-16) 

(4-17) 

(4-18) 

(4-19) 

re- 

( 4-20) 

linear 


equations  can  be  described  in  matrix  terms  by 


corresponding  Co 


X  a 


(4-22) 


where 


a  = 


kx^_1(l  ♦  x  )q  ♦  x^q(  1  +  x,)*1'1  r  n 

- 2  C" .  . .  "  2L  X2(l  +  X2]  J  (4“23) 


We  again  derive  an  expression  for  A  in  terns  of  x^  by  using  deter¬ 
minants*  and  solve  for  the  value  of  Xj  that  results  in  A  =0.  This  equa¬ 
tion  is 
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I  xf(l  +  X  )q  x*(l  +  x  )q  | 

■  K  -  *»>  *  1  ct  ‘  -  I  4 


( 4-24) 
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,  ,  (x2-x2)x*(l+x.)q  x2x*(H-x.)q  x«x![(  1+x,  )q  j 

♦  -  1  V  -1\  3  -  :LVX  I 
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where  a  is  as  defined  in  eq  (4-23). 

Equation  (4-24)  is  solved  numerically*  and  the  resulting  value  of  x2 

is  used  to  determine  the  corresponding  maximum  realizable  value  of  D. 
Again  using 


D  = 


Num( x- ) 
Den(x2) 


(4-25) 


we  have*  for  the  distributed  case 


Num(x, 


]  ^  *  y  j 
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L 


and 


Den(x2)  =  (xj-x3)x2  ♦  2(xj-Xj)x2  ♦  (xj-3xjX3-3XjX3~x3)x2  (4-26) 

,  33  .  .32  2  3. 

+  2(XjX3-XjX3)x2  ♦  (XjXj-XjXj) 
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After  substituting  the  known  values  into  equations  (4-26)  and  (4-27) 
the  resulting  maximum  realizable  value  for  D  can  be  calculated  via 
eq  (4-25). 

This  procedure  has  also  been  coded  into  a  set  of  computer  subroutines 
(listed  in  Appendix  A),  and  incorporated  into  the  ADIST  computer  program 
for  distributed  network  synthesis  where  the  maximum  allowable  DPC  value 
for  3rd  order  networks  is  calculated  and  presented  to  the  operator  at  the 

appropriate  time. 


Chapter  5 


Restricted  Topology  Matching  letworks 
for  Bose  of  Fabricatioo 


When  experimenting  with  different  structures  and  doping  profiles  to 
optimize  FET's  intended  for  monolithic  applications*  it  is  convenient  to 
have  impedance  matching  topologies  available  that  require  no  capacitors 
and  only  one  level  of  metallization*  These  circuits  would  make  it  possi¬ 
ble  to  quickly  and  easily  fabricate  single-stage  monolithic  amplifiers  to 
realistically  evaluate  FET  performance*  Since  these  test  amplifiers  will 
require  off-chip  biasing  (usually  through  bias  tees)*  their  use  will  be 
largely  limited  to  the  laboratory  —  but  there  they  can  be  very  useful  for 
demonstrating  the  capabilities  of  a  particular  FET  design  in  a  realistic 
amplifier  environment. 


i*l*  Allowable  Topologies 


The  circuit  topologies  which  meet  the  restrictions  described  above 
are  composed  only  of  shunt  open-circuited  stubs  and  cascaded  transmission 
line  sections*  Matching  networks  composed  only  of  these  two  circuit  ele¬ 
ments  require  only  one  level  of  metallization*  no  RF  shorts  to  the  chip 
ground  plane*  and  no  capacitors*  These  traits  make  this  class  of  circuit 
very  desirable  for  some  applications*  but  there  are  also  disadvantages* 
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In  addition  to  the  problem  of  bias  introduction,  these  circuits  allow 
only  indirect  control  of  the  impedance  transformation  ratio  provided  by 
the  network.  Since  shunt  short-circuited  stubs  are  not  allowed  in  these 
topologies*  the  high-pass  Kuroda  Transformation  (§  3.8)  cannot  be  used  to 
adjust  impedance  transformation  ratio.  As  a  result*  the  task  of  simul¬ 
taneously  providing  a  specific  impedance  transformation  while  exactly 
absorbing  the  FET  parasitic  reactance  is  much  more  difficult  with  these 
networks.  Indeed,  in  many  cases  no.  directly  synthesized  network  will  be 
found  that  can  simultaneously  absorb  the  FET  parasitic  and  match 
impedances.  In  these  cases*  however*  using  a  cascade  of  two  networks  will 
usually  solve  the  problem*  though  sometimes  requiring  computer  optimiza¬ 
tion  for  broadband  applications. 

Input  networks. 

Several  input  network  topologies  that  conform  to  these  topology  res¬ 
trictions  are  shown  in  Figure  5.1.  The  3rd  order  circuit  shown  at  Figure 
5.1(a)  may  be  synthesized  via  the  ADIST  computer  program  for  broadband 
applications,  or  by  Smith  Chart  for  narrowband  applications.  The  4th  and 
5th  order  topologies  shown  at  Figure  5.1(b-d)  can  be  directly  synthesized 
via  ADIST  for  either  wide  or  narrowband  applications*  while  the  circuit  at 
Figure  5.1(e)  is  formed  by  cascading  4^  and  2nd  order  sections. 


The  cascading  procedure  here  is  similar  to  that  described  in  |  2.9  and 
8  3.6;  this  technique  is  even  more  important  when  no  directly  syntheaiaed 
network  can  be  found  that  simultaneously  absorbs  the  device  parasitic  and 
produces  the  necessary  impedance  transformation. 

1.1.2.  Output  networks 

Usable  output  topologies  composed  only  of  shunt  open  stubs  and  cas¬ 
cade  lines  are  shown  in  Figure  5.2.  The  topology  shown  at  Figure  5.2(a) 
may  be  synthesized  by  Smith  Chart  for  narrowband  applications*  and  those 
at  (b)  and  (c)  may  be  synthesized  by  the  ADIST  program  for  either  narrow 
or  wide  bandwidths.  The  topology  shown  in  Figure  5.2(d)  is  formed  by  cas¬ 
cading. 

1.2.  Single-Frequency  Oeaign  Method 

The  simple  networks  of  Figures  5.1(a)  and  5.2(a)  may  be  designed  for 
a  match  at  a  single-frequency  by  using  the  Smith  Chart.  Figure  5.3  graph¬ 
ically  illustrates  the  procedure. 

Starting  from  the  normalized  FET  input  or  output  impedance  (point  A)* 
a  cascade  transmission  line  section  is  used  to  proceed  to  point  B.  The 
length  of  this  line  is  chosen  so  that  after  reflecting  point  B  into  the 
admittance  plane  (at  point  C),  C  lines  on  the  g=l  circle.  The  match  is 
then  completed  by  moving  from  C  to  D  by  adding  a  shunt  open-circuited 


£•1.  Broadband  Design  Method 

£.1.1.  Direct  Synthesis 

It  is  always  best  to  use  a  directly  synthesized  broadband  network 
(rather  than  a  cascaded  design)  when  a  suitable  combination  of  topology 
and  gain  function  specifications  can  be  found  to  simultaneously  absorb  the 
FET's  parasitic  capacitance  and  provide  the  necessary  impedance  transfor¬ 
mation.  With  the  restricted  topology  networks  of  this  section  it  is  espe¬ 
cially  important  that  all  possible  topologies,  bandwidths.  ripple  specifi¬ 
cations.  and  Sjj(p)  zero  distributions  be  explored  when  seeking  a  solu¬ 
tion.  since  Kuroda  Transformations  may  not  be  used. 

Due  to  this  lack  of  flexibility,  in  many  cases  it  will  not  be  possi¬ 
ble  to  find  a  suitable  directly  synthesized  network.  In  these  cases,  a 
cascaded  design  may  be  used. 

1.1.2.  Cascaded  Synthesis 

£.1.1.1.  Flat  Match 

The  cascading  procedure  when  an  impedance  match  with  no  gain  slope 
compensation  is  desired  is  as  follows: 

(1)  Synthesize  a  zero  gain  slope  3rd  or  4th  order  restricted 
topology  network  that  meets  the  gain  and  bandwidth  requirements 
while  exactly  absorbing  the  FET  parasitic  input  (or  output)  capaci¬ 
tance  —  but#without  regard  to  the  required  network  terminating 
impedance  (R  ). 
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(2)  Synthesize  a  zero  gain  slope  2nd  order  restricted  topology 
netvork  that*  over  thg  specified  band*  produces  an  impedance 
transformation  from  R  to  the  desired  level  at  the  source  (or  load) 
for  the  amplifier. 

(3)  Cascade  these  two  circuits*  analyze  the  resulting  gain 
response*  and  apply  computer  optimization*  if  required*  to  yield  a 
satisfactory  gain  response. 


This  procedure  yields  circuits  like  those  shown  in  Figures  5.1(e)  and 
5.2(d).  Very  satisfactory  results  have  been  obtained  with  this  technique 
up  to  octave  bandwidth. 


1.1.2,*!.  fiAin.  Sloped  Match 


The  cascading  procedure  for  an  impedance  match  with  gain  slope  com¬ 
pensation  is  slightly  different  than  that  described  above*  and  it  is  not 
possible  to  exactly  synthesize  a  specified  value  of  gain  slope  compensa¬ 
tion.  Due  to  mismatch  effects  between  the  two  cascaded  networks  used*  the 
final  slope  of  the  cascade  is  impossible  to  predict  exactly.  An  approxi¬ 
mate  procedure  is  therefore  used*  and  one  or  two  iterations  of  the  pro¬ 
cedure  may  be  required  to  obtain  the  desired  slope  in  the  final  network. 
The  procedure  is  as  follows: 


(1)  Synthesize  a  3rd  or  4th  order  restricted  topology  network 
that  provides  one-half  the  desired  gain  slope  compensation  over  the 
required  band  while  exactly  absorbing  the  FET  parasitic  capacitance 
—  again  without  regard  to  the  required  network  terminating 
impedance  (R  ) . 

(2)  Synthesize  a  zero  gain  slope  2nd  order  restricted  topology 
network  that  over  twice  the  specified  hand  produces  an  impedance 
transformation  from  R*  to  the  desired  level.  (This  network  is 
designed  over  twice  the  specified  band*  with  its  center  frequency  at 
the  upper  edge  of  the  specified  passband.  In  this  way*  an  exact 
impedance  match  is  provided  at  the  upper  passband  edge  with  increas¬ 
ing  insertion  loss  at  lower  frequencies.) 
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(3)  Cascade  these  two  circuits  and  apply  computer  optimization 
to  obtain  a  satisfactory  gain  slope* 

(4)  If  necessary,  increase  or  decrease  the  starting  value  of 
gain  slope  used  in  step  (1)  and  repeat  the  procedure  to  obtain  a 
more  satisfactory  result. 


i*A*  Restricted  Topology  Design  Examples 


Narrow  and  broadband  monolithic  FET  power  amplifiers  were  designed  by 
using  the  above  procedures,  and  the  resulting  amplifiers  were  fabricated 
on  GaAs  and  tested  by  another  student.  Details  of  the  device  characteri¬ 
zation.  GaAs  fabrication,  and  RF  measurement  procedures  used  can  be  found 
in  [12]. 


1.4.1*  Octave  Rand.  Amplifier 

Figure  3.4  shows  the  6-12  GHz  monolithic  power  amplifier  designed  by 
using  these  techniques.  The  power  FET  was  a  600  pm  gate  width  design  with 
a  0.8  pm  gate  length  (additional  details  on  the  device  can  also  be  found 
in  [12]).  The  figure  shows  the  FET  modelled  unilaterally  for  design  of 
the  matching  networks.  The  input  model  was  obtained  by  optimization  to 
the  measured  small-signal  8^  characteristics  of  the  device,  while  the 

output  model  was  optimized  to  model  the  impedance  that  maximized  output 
power  from  the  device,  as  measured  under  large-signal  conditions. 

The  input  network  shown  in  Figure  3.4  was  formed  by  producing  a 
single-frequency  match  at  12  GHz.  the  upper  band  edge.  Analysis  showed 
that  this  circuit  yielded  a  smoothly  decreasing  gain  response  with  7.3  dB 
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Figure  5.4.  Restricted  Topology  Monolithic  6-12  GHz  Power  Amplifier 


of  insertion  loss  st  the  lower  band  edge.  6  GHz  --  very  satisfactory  for 
gain  slope  compensation  of  this  device.  No  computer  optimization  of  this 
network  was  required . 

The  output  network  was  formed  by  cascading  41*1  and  2n<*  order  sec¬ 
tions;  it  was  designed  for  zero  gain  slope.  (No  directly  synthesized 
network  meeting  the  topology  restrictions  of  this  section  could  be  found 
to  exactly  match  the  device  to  a  50  Ohm  load  over  the  6-12  GHz  band.) 

The  FET's  parasitic  output  capacitance  was  absorbed  by  a  4th  order  network 
(Z^,  Z2.  and  Z^  in  conjunction  with  the  FET's  capacitance);  this  network 
required  a  5.2  Ohm  terminating  impedance. 

A  2n<*  order  network  consisting  of  the  line  section  Z^  and  an  open 
shunt  stub  (not  shown)  was  then  synthesized  to  accomplish  the  5.2  to  50 
Ohm  impedance  transformation  over  6-12  GHz.  The  two  networks  were  cas¬ 
caded  and  computer  optimization  was  applied  to  minimize  insertion  loss 
across  the  band.  This  optimization  reduced  the  open  shunt  stub  to  an 
insignificant  line  length;  the  stub  was  then  removed  from  design  and 
a  final  optimization  performed. 

The  predicted  and  measured  performance  of  the  final  amplifier  circuit 
is  also  shown  in  Figure  5.4.  The  amplifier  achieved  5  i  0.5  dB  power  gain 
over  6-12  GHz  when  operating  at  the  design  drive  level  of  +17  dBm. 

Figure  5.5  shows  the  gain  performance  of  the  cascaded  output  network 
before  and  after  computer  optimization. 


INSERTION  LOSS  (dB) 
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Figure  5.6  is  a  scanning  electron  microscope  (SEM)  photograph  of  the 
6  -  12  GHz  monolithic  amplifier  chip  prior  to  scribing.  (The  distortion 
apparent  in  the  image  is  due  to  irregularities  in  the  electron  beam  scan 
at  the  extreme  angle  required  to  image  this  particular  chip  on  the  wafer.) 


i.4t*2.  JJ2X  Band  Amplifier 

Figure  5.7  shows  the  8.5  -  9.5  GHz  narrowband  monolithic  power 
amplifier  designed  by  these  techniques.  The  FET  was  the  same  as  that  used 
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Figure  5.5.  Wideband  Output  Network  Performance 


Figure  5.6.  SEM  Photograph  of  Octave  Band  Power  Amplifier 


in  the  octave  band  design •  above. 

The  input  network  shown  in  Figure  5*7  was  formed  by  producing  a 
single-frequency  match  at  9.5  GHz*  the  upper  band  edge.  Analysis  showed 
that  this  networks  too*  provided  adequate  gain  slope  compensation  for  the 
FET.  Ho  computer  optimization  was  required  for  this  network. 

The  output  network  was  formed  by  producing  a  single-frequency  match 
at  9  GHz*  the  band  center.  This  network  exhibited  0.1  dB  insertion  loss 
at  the  band  edges  and  0  insertion  loss  at  9  GHz;  no  computer  optimization 
was  required. 
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Figure  S.7.  Restricted  Topology  Monolithic  10Z  Band  Power  Amplifier 
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Figure  5*7  also  shows  the  calculated  gain  response  of  this  amplifier 
based  on  the  measured  large-signal  power  contours  for  the  FET.  This 
amplifier  achieved  a  (predicted)  7.5  i  0.5  dB  power  gain  at  its  design 
input  drive  level  of  +13. 5  dBm.  Fabrication  of  this  amplifier  was  not 
successful  due  to  processing,  so  measured  results  for  this  design  are  not 
available. 

Figure  5.8  is  a  scanning  electron  microscope  photograph  of  the 
8.5  -  9.5  GHz  amplifier  as  fabricated. 


Figure  5.8.  SEM  Photograph  of  Monolithic  101  Band  Power  Amplifier 


Chapter  6 


Cornell 


The  computer-aided  procedures  for  systematic  design  of  lumped  element 
and  distributed  element  amplifier  matching  networks  (limited  to  even 
order)  have  been  extended  to  allow  inclusion  of  odd-order  networks.  The 
synthesis  method  used  for  these  odd-order  networks  introduces  a  new  syn¬ 
thesis  parameter*  the  denominator  polynomial  constant  (DPC)»  which  pro¬ 
vides  an  added  degree  of  freedom  for  the  circuit  designer*  making  it 
easier  to  meet  specified  circuit  performance  goals  with  very  simple  match¬ 
ing  networks  (such  as  3r<*  order).  The  final  synthesis  algorithm  allows 

specification  of  either  even-  or  odd-order  network  topologies  for  practi¬ 
cal  computer  implementation. 

The  gain-bandwidth  performance  of  general  odd-order  networks  wa6  dis¬ 
cussed*  and  specific  gain-bandwidth  results  were  presented  for  several 
different  circuit  topologies  as  examples.  An  algorithm  for  determining 
the  DPC  value  which  maximizes  the  gain-bandwidth  product  of  an  odd-order 
network  was  presented*  and  practical  circuit  topologies  for  monolithic 
input*  output*  and  interstage  networks  were  discussed  along  with  their 
design  procedures  for  both  lumped  and  distributed  element  circuits. 

The  very  significant  decrease  in  monolithic  matching  network  complex¬ 
ity  made  possible  by  the  use  of  odd-order  networks  was  demonstrated  by 
examples  where  3  element  networks  were  designed  to  replace  much  larger  5 
element  matching  networks*  usually  without  significant  performance  degra- 
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datioo.  Both  low-noise  and  power  amplifier  designs  were  illustrated. 

This  improvement  was  made  possible  by  choosing  DPC  values  for  the  3rd 
order  networks  so  that  the  extra  inductor  normally  required  to  provide  a 
specific  impedance  transformation  ratio  by  a  Norton  Transformation  in  the 
network  could  be  eliminated. 

An  analytical  technique  was  developed  for  determining  the  maximum 
realizable  value  of  the  denominator  polynomial  constant  in  3rd  order  gain 

functions  for  both  lumped  and  distributed  elements,  regardless  of  ripple, 
bandwidth,  minimum  insertion  loss,  or  gain  slope  specifications.  This 
information  is  useful  to  the  designer  at  the  time  of  synthesis,  and  com¬ 
puter  implementations  of  the  procedures  were  presented  that  have  been 
incorporated  into  network  synthesis  programs. 

Chapter  5  presented  design  techniques  and  practical  topologies  for 
matching  networks  using  only  cascaded  lines  and  shunt  open-circuited  stubs 
for  extremely  simple  monolithic  realizations.  A  significant  limitation 
with  these  circuits  is  the  very  limited  control  available  over  impedance 
transformation  ratio;  this  problem  was  resolved  by  developing  a  cascaded 
synthesis  technique  which  includes  final  computer  optimization. 
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Appendix  A 


Source  Code  Listings 


FORTRAN  source  code  listings  for  the  following  subroutines  are 
included  in  this  appendix: 

A.l  ODDSYN  (Lumped  element  gain  function  synthesis) 

A. 2  ODDIST  (Distributed  element  gain  function  synthesis) 

A. 3  DLIMTL  (Lumped  element  max  DPC  calculation) 

A. 4  DLIMITD  (Distributed  element  max  DPC  calculation) 
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ODDSTB 


A..1  QDDSYN 

SUBROUTINE  ODDSYN  (  N,  XX,  RIP,  GS.  BWR,  K.  10,  DPC  ) 

C  J.  T.  DIJAK  16  MAR  1983 

DOUBLE  PRECISION  BWR,  RIP,  GS,  X(7),  GSP(7).  A(7.7).  DPC. 

1  F(7).  C(7) ,  XX(22) ,  WO,  Wl,  Z,  STEP,  XINC,  GMAX 

INTEGER  I,  J,  K.  N,  IPVT(7),  FLAG.  P(6,6).  PASS.  10 
LOGICAL  ODD,  I FLAG 

DATA  (P( 1 ,J) ,J=1 .6)  /  O.O.O.O.O.O  / 

DATA  (P(2 ,J) ,J=1 .6)  /  1.2,0 ,0,0.0  / 

DATA  (P(3,J).J=1.6)  /  2. 1,3 ,0,0.0  / 

DATA  (P(4,J) ,J=1 .6)  /  1.4, 2. 3.0,0  / 

DATA  (P(5.J).J=1.6)  /  2, 5,1.3, 4,0  / 

DATA  (P(6 <J) ,J=1 ,6)  /  1.6, 3. 4, 2, 5  / 

IF  (N  .GT.  6)  N  =  6 
IF  (N  .LT.  2)  N  =  2 

C  SET  ODD  ORDER  FLAG  IF  N  ODD. 

ODD  =  .TRUE. 

IF  (N  .EQ.  (2*((N+l)/2)))  ODD  =  .FALSE. 

DO  5  I  =  1.  N+l 
C(N)  =  1.0D0 
5  CONTINUE 

C  SET  UP  VECTOR  OF  CRITICAL  FREQUENCIES. 

X(l)  =  1.0D0  /  (BWR  *  BWR) 

IF  (ODD)  GO  TO  15 
X(  N+l  )  =  1.0D0 
N2  =  N 

XINC  =  (1.0D0  -  X(l))  /  DFLOAT(N) 

DO  10  I  =  2,  N2 

X(I)  =  X(I-l)  +  XINC 
10  CONTINUE 

GO  TO  25 

C  ODD-  ORDER  CRITICAL  FREQUENCIES. 

15  X(N)  =  1.0D0 

N2  =  N  -  1 

XINC  =  (1.0D0  -  X(l))  /  DFLOAT(N-l) 

DO  20  I  =  2,  N2 

X(I)  =  X(I-l)  +  XINC 
20  CONTINUE 

C  SET  (N+1)ST  POLYNOMIAL  COEFFICIENT  IF  ODD. 

C(N+1)  =  DPC 

IF  (  N  .LT.  4  )  GO  TO  30 

X(2)  =  X(2)  -  XINC  /  2.0D0 
X(N2)  =  X(N2)  +  XINC  /  2.0D0 
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30  N3  =  N  +  1 

IF  (  ODD  )  N3  =  N 

IF  (10  .GT.  0)  WRITE  (6.920)  (  X(I).  I  =  1.  N3  ) 

920  FORMAT  (  /  '  INITIAL  CRITICAL  FREQUENCIES  /  7F10.5  ) 

PASS  =  1 

C  SET  UP  THE  A  MATRIX  IN  SCRAMBLED  FORM. 

35  DO  40  I  =  1.  N 

DO  40  J  =  1.  N 

A(I.J)  =  X(I)  **  P(N,J) 

40  CONTINUE 

IF  (ODD)  GO  TO  65 
DO  45  I  =  1.  N3 
A(I.N3)  =  l.ODO 
A(N3.I)  =  l.ODO 
45  CONTINUE 

C  COMPUTE  SPECIFIED  GAIN  AT  CRITICAL  FREQUENCIES. 

65  IF  (  GS  ,NE.  0.0  )  GO  TO  75 

DO  70  I  =  1.  N3 

GSP(I)  =  DEXP(  -RIP  *  0.23025851D0  )  * 

2  DEXP(  RIP  *  0.23025851D0  *  (-1)**I) 

70  CONTINUE 

GO  TO  85 

75  GSP(N3)  =  DEXP(  -2.0D0  *  RIP  *  0.23025851D0  ) 

N2  =  N3  -  1 

DO  80  I  =  1.N2 

GSP(I)  =  GREF(  X(I),  GS  ) 

IF  ((I.EQ.l) .0R.(I.EQ,3) .0R.(I.EQ.5))  GSP(I)  =  GSP(I)*GSP(N3) 
80  CONTINUE 

C  COMPUTE  RIGHT  HAND  SIDE  OF  MATRIX  EQUATION. 

85  Z  =  O.ODO 

IF  (  ODD  )  Z  =  C(N+1) 

DO  90  I  =  1.N3 

IF  (K  .EQ.  0)  F( I)  =  l.ODO  /  GSP(I)  -  Z 
IF  (K  .NE.  0)  F(I)  =  (  X( I ) **K  )  /  GSP(I)  -  Z 
90  CONTINUE 

C  DO  LU  FACTORIZATION  OF  MATRIX  A. 

CALL  DGEFA  (  A,  7.  N3,  IPVT,  FLAG  ) 

IF  (FLAG  .EQ.  0)  GO  TO  110 
WRITE  (6,930) 

930  FORMAT  (  /  '  LU  FACTORIZATION  FAILURE.'  /  ) 

RETURN 

C  SOLVE  FOR  VECTOR  OF  COEFFICIENTS. 

110  CALL  DGESL  (  A,  7,  N3,  IPVT,  F,  0  ) 

C  UNSCRAMBLE  POLYNOMIAL  COEFFICIENTS  INTO  DESCENDING  ORDER. 
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DO  115  I  =  1.  N 

C(  N+1-P(N.I)  )  =  F(l) 

115  CONTINUE 

IF  (.NOT.  ODD)  C(N+1)  =  F(N+1) 

C  FIND  TRUE  CRITICAL  FREQUENCIES. 

DO  150  I  =  2.  N2 

C  RESET  FLAG  SO  WE  CAN  DETERMINE  IF  ANY  INTERATIONS  WERE  NEEDED. 

I FLAG  =  .FALSE. 

IF  (PASS  .LE.  2)  WO  =  X(I-l) 

IF  (PASS  .GT.  2)  WO  =  X(I) 

DO  1 40  J  =  1.  10 

120  IF  ((I.EQ.2)  .OR.  (I.EQ.4)  .OR.  (I.EQ.6))  GO  TO  130 

125  IF  (GPP(WO.C.K.N)  .GT.  O.ODO)  GO  TO  135 

WO  =  1.02D0  *  WO 
GO  TO  125 

130  IF  (GPP(WO.C.K.N)  .LT.  O.ODO)  GO  TO  135 

WO  =  1.02D0  *  WO 
GO  TO  130 

135  STEP  =  (GREFP(WO.GS)  -  GP(WO.C.K.N))  / 

2  (GREFPP(WO.GS)  -  GPP(WO.C.K.N) ) 

IF  (STEP  .LT.  -0.03D0)  STEP  =  -0 .03D0 
IF  (STEP  .GT.  +0.03D0)  STEP  =  +0.03D0 
W1  =  WO  -  STEP 

Z  =  GREFP (  Wl.GS  )  -  GP(  Wl.C.K.N  ) 

IF  (DABS(  Z  )  .LT.  0.00005D0)  GO  TO  145 
WO  =  W1 

I  FLAG  =  .TRUE. 

140  CONTINUE 

145  X( I)  =  W1 

150  CONTINUE 

C  HERE  ON  TERMINATION  OF  LOOP  TO  FIND  TRUE  CRITICAL  FREQUENCIES. 

IF  (10  .GT.  0)  WRITE  (6,940)  (X(I),I  =1,  N3) 

940  FORMAT  (  /  '  RE-COMPUTED  CRITICAL  FREQUENCIES  :*  /  7F10.5  ) 

PASS  =  PASS  +  1 

IF  (.NOT.  I FLAG)  GO  TO  160 

IF  (PASS  .LT.  10)  GO  TO  35 

C  JUMP  HERE  AFTER  COMPLETING  FINAL  PASS  THROUGH  SYNTHESIS. 

160  N2  =  N  +  1 

C  LOAD  DENOMINATOR  COEFFICIENTS  INTO  OUTPUT  VECTOR. 

DO  165  1=1.  N2 

XX(I)  =  C(N2  +1-1) 

165  CONTINUE 

C  LOAD  CRITICAL  FREQUENCIES  INTO  OUTPUT  VECTOR. 

DO  170  1=1.  N3 

XX(N2  +  I)  =  X(I) 

170  CONTINUE 

IF  (10  .GT.  0  )  WRITE  (6,980)  (XX(I),  1=1,  N+l  ) 

980  FORMAT  (/  '  DENOMINATOR  COEFFICIENTS:  (Constant  first)'  /  7F11.6  ) 
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IF  (10  .LT.  0)  RETURN 

C  FIND  MAX  VALUE  OF  GAIN  FUNCTION,  SO  WE  CAN  NORMALIZE  TO  1.0. 
CALL  MXGAIN(  GMAX.X.C.K.N  ) 

C  NORMALIZE  GAIN  FUNCTION,  AND  COMPUTE  RESPONSE. 

WRITE  (6,950) 

950  FORMAT  (  /  T17 .  'GAIN*  /  T7 .  'FREQ',  T17.  '(DB)'  ) 

XINC  =  (1.0D0  -  X(l))  /  10.0D0 
DO  180  I  =  1,  11 
C  IN  BAND  RESPONSE. 

W1  =  X(l)  +  FLOAT(I-l)  *  XINC 

G1  =  10 .0D0*DLOG10(  GAIN2(  Wl.C.K.N  )  /  GMAX  ) 

W2  =  DSQRT(  W1  ) 

WRITE  (6,960)  W2,  G1 

960  FORMAT  (  T4.  F7.3.  F10.3  ) 

180  CONTINUE 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GAIN 2  (  W,  C.  K,  N  ) 

DOUBLE  PRECISION  W.  C(7),  Z 
INTEGER  N,  K 
Z  =  W  *  C(l) 

DO  10  I  =  2.  N 

Z  =  W  *  (Z  +  C(D) 

10  CONTINUE 

Z  =  Z  +  C(N+1) 

IF  (K  .EQ.  0)  GAIN 2  =  1 .0D0  /  Z 
IF  (K  .NE.  0)  GAIN2  =  (W**K)  /  Z 
RETURN 
END 

SUBROUTINE  MXGAIN (  GMAX,  X,  C,  K,  N  ) 

C  FIND  MAXIMUM  GAIN  POINT  OVER  X(l)  TO  2  -  X(l)  FREQS. 

DOUBLE  PRECISION  X(7),  C(7),  W,  GMAX,  XINC,  G 
INTEGER  N,  K 

XINC  =  (1.0D0  -  X(l))  /  40.0D0 
W  =  X(l) 

GMAX  =  -25.0D0 
DO  10  I  =  1,  80 
W  =  W  +  XINC 
G  =  GAIN 2(  W,  C,  K,  N  ) 

IF  (G  .GT.  GMAX  )  GMAX  =  G 
10  CONTINUE 

RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GP  (  W,  C,  K.  N  ) 

C  COMPUTE  VALUE  OF  FIRST  DERIVATIVE  OF  GAIN  FUNCTION  AT  W. 

DOUBLE  PRECISION  W,  C(7).  W1 ,  W2 
INTEGER  K .  N 
W1  =  W  ♦  0.0005D0 
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W2  =  W  -  0.0005D0 

GP  =  GAIN2(W1,C.K.N)  /  0.001D0  -  GAIN2(W2.C.K.N)  /  0.001D0 
RETURN 
END 


DOUBLE  PRECISION  FUNCTION  GPP  (  W,  C.  K.  N  ) 

COMPUTE  VALUE  OF  SECOND  DERIVATIVE  OF  GAIN  FUNCTION  AT  W. 
DOUBLE  PRECISION  W,  C(7),  W1 ,  W2 
INTEGER  K.  N 
W1  =  W  +  0.0005D0 
W2  =  W  -  0.0005D0 

GPP  =  GP(W1,C,K,N)/0.001D0  -  GP(W2.C.K.N)/0.001DO 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GREF  (  W,  GS  ) 

COMPUTE  REFERENCE  GAIN  AT  W. 

DOUBLE  PRECISION  W,  GS 

GREF  =  DEXP(  0.38245111D0  *  GS  *  DLOG10(  W  )  ) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  GREFP  (  W.  GS  ) 
COMPUTE  FIRST  DERIVATIVE  OF  GREF(X)  AT  W. 
DOUBLE  PRECISION  W,  GS 

GREFP  =  0.1660964D0  *  GS  *  GREF(  W,  GS  )  /  W 
RETURN 
END 


DOUBLE  PRECISION  FUNCTION  GREFPP  (  W,  GS  ) 

C  COMPUTE  SECOND  DERIVATIVE  OF  GREF(X)  AT  W. 

DOUBLE  PRECISION  W.  GS 

GREFPP  =  (GS  *  0. 166096 4D0  /  W)*(  GREFP(W.GS)  -  GREF(W.GS)/W  ) 
RETURN 
END 
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A.2  QHHLSlI 

SUBROUTINE  ODDIST  (  N.  XX.  RIP.  GS.  NPL .  K,  10.  DPC  ) 

C  THIS  ROUTINE  IS  COMPATIBLE  WITH  THE  DIST  PROGRAM  AND 
C  TAKES  THE  PLACE  OF  THE  MODNWT  ROUTINE.  ODDIST  SYNTHESIZES  EVEN  OR 
C  ODD  ORDER  GAIN  FUNCTIONS  OF  ORDER  2  THROUGH  6  BY  ITERATIVELY  SOLVING 
C  A  SYSTEM  OF  LINEAR  EQUATIONS  UNTIL  BOTH  THE  GAIN  PERFORMANCE  AND 
C  THE  CRITICAL  FREQUENCIES  OF  THE  SYNTHESIZED  FUNCTION  ARE  KNOWN. 

C  ODDIST  DOES  FOR  DISTRIBUTED  CIRCUITS  WHAT  ODDSYN  DOES  FOR  LUMPED. 

C  J.  T.  DIJAK  16  MAY  1983 

DOUBLE  PRECISION  RIP.  GS ,  X(7).  GSP(7).  A(7.7).  DPC. 

1  F(7) .  C(7).  XX( 22) .  WO.  W1 .  Z,  STEP.  XINC.  GMAX 

INTEGER  I,  J,  K.  N.  IPVT(7),  FLAG,  P(6.6),  PASS,  10.  NPL 
LOGICAL  ODD.  I FLAG 

DATA  (P(l ,J) ,J=1 ,6)  /  0. 0,0.0. 0.0  / 

DATA  (P(2.J) ,J=1 ,6)  /  1 ,2.0.0 .0,0  / 

DATA  (P(3.J) ,J=1 ,6)  /  2. 1,3. 0.0.0  / 

DATA  (P(4,J) ,J=1 ,6)  /  1.4.2.3.0.0  / 

DATA  (P(5.J) ,J=1 .6)  /  2, 5. 1.3. 4.0  / 

DATA  (P(6 ,J) ,J=1 .6)  /  1.6,3. 4.2. 5  / 

IF  (N  .GT.  6)  N  =  6 
IF  (N  .LT.  2)  N  =  2 

C  SET  ODD  ORDER  FLAG  IF  N  ODD. 

ODD  =  .TRUE. 

IF  (N  .EQ.  ( 2*( (N+l)/2) ) )  ODD  =  .FALSE. 

DO  5  1=1.  N+l 

C(N)  =  1.0D0 
5  CONTINUE 

C  SET  UP  VECTOR  OF  CRITICAL  FREQUENCIES. 

X(l)  =  XX(N+2) 

IF  (ODD)  GO  TO  15 
X(  N+l  )  =  XX( 2*N+2) 

N2  =  N 

XINC  =  (X(N+1 )  -  X(l))  /  DFLOAT(N) 

DO  10  I  =  2.  N2 

X(I)  =  X(I-l)  +  XINC 
10  CONTINUE 

GO  TO  25 

C  ODD-  ORDER  CRITICAL  FREQUENCIES. 

15  X(N)  =  XX( 2*N+2) 

N2  =  N  -  1 

XINC  =  (X(N)  -  X(l))  /  DFLOAT(N-l) 

DO  20  I  =  2.  N2 

X(I)  =  X(I-l)  +  XINC 
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20  CONTINUE 

C  SET  (N+1)ST  POLYNOMIAL  COEFFICIENT  IF  ODD. 

C(N+1 )  =  DPC 

25  IF  (  N  .LT.  4  )  GO  TO  30 

X(2)  =  X(2)  -  XINC  /  2.0D0 
X(N2)  =  X(N2)  +  XINC  /  2.0D0 

30  N3  =  N  +  1 

IF  (  ODD  )  N3  =  N 

IF  CIO  .GT.  0)  WRITE  (6.920)  (  X(I).  I  =  1 .  N3  ) 

920  FORMAT  (  /  '  INITIAL  CRITICAL  FREQUENCIES  :»  /  7F10.5  ) 

PASS  =  1 

C  SET  UP  THE  A  MATRIX  IN  SCRAMBLED  FORM. 

35  DO  40  I  =  1.  N 

DO  40  J  =  1.  N 

A(I.J)  =  X(I)  **  P(N.J) 

40  CONTINUE 

IF  (ODD)  GO  TO  65 
DO  45  I  =  1.  N 
A(I.N3)  =  1.0D0 
A(N3.I)  =  X(N3)**P(N.I) 

45  CONTINUE 

A(N3.N3)  =  1.0D0 

C  COMPUTE  SPECIFIED  GAIN  AT  CRITICAL  FREQUENCIES. 

65  IF  (  GS  .NE.  0.0  )  GO  TO  75 
DO  70  I  =  1.  N3 

GSP(I)  =  DEXP(  -RIP  *  0.23025851D0  )  * 

2  DEXP(  RIP  *  0.23025851D0  *  (-1)**I) 

70  CONTINUE 

GO  TO  85 

75  Z  =  DEXP(  -2.0D0  *  RIP  *  0.23025851D0  ) 

N2  =  N3  -  1 

DO  80  I  =  1.N3 

GSP(I)  =  GREF(  X(I).  GS  ) 

IF  ((I.EQ.l) .OR.(I.EQ.3)  .OR.(I.EQ.S))  GSP(I)  =GSP(I)*Z 
80  CONTINUE 

C  COMPUTE  RIGHT  HAND  SIDE  OF  MATRIX  EQUATION. 

85  Z  =  0.0D0 

IF  (  ODD  )  Z  =  C(N+1) 

DO  90  I  =  1.N3 

IF  (K  .EQ.  0)  GO  TO  88 

IF  (NPL  .EQ.  0)  F(I)  =  (X(I)**K)/GSP(I)  -  Z 
IF  (NPL  .NE.  0)  F( I ) =( (X( I )**K)*(X(I)+l .DO)**NPL)/GSP(I)-Z 
GO  TO  90 

IF  (NPL  .EQ.  0)  F(I)  =  1 .DO  /  GSP(I)  -  Z 

IF  (NPL  .NE.  0)  F(I)  =  (X(I)  +  1 .D0)**NPL  /  GSP(I)  -  Z 
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C  DO  LU  FACTORIZATION  OF  MATRIX  A. 

CALL  DGEFA  (  A.  7.  N3,  IPVT.  FLAG  ) 

IF  (FLAG  .EQ.  0)  GO  TO  110 
WRITE  (6.930) 

930  FORMAT  (  /  •  LU  FACTORIZATION  FAILURE.'  /  ) 

RETURN 

C  SOLVE  FOR  VECTOR  OF  COEFFICIENTS. 

110  CALL  DGESL  (  A.  7.  N3 .  IPVT.  F.  0  ) 

C  UNSCRAMBLE  POLYNOMIAL  COEFFICIENTS  INTO  DESCENDING  ORDER. 

DO  115  I  =  1.  N 

C(  N+1-P(N,I)  )  =  F(l) 

115  CONTINUE 

IF  (.NOT.  ODD)  C(N+1 )  =  F(N+1) 

C  FIND  TRUE  CRITICAL  FREQUENCIES. 

DO  150  1=2.  N2 

C  RESET  FLAG  SO  WE  CAN  DETERMINE  IF  ANY  ITERATIONS  WERE  NEEDED. 

I FLAG  =  .FALSE. 

IF  (PASS  .LE.  2)  WO  =  X(I-l) 

IF  (PASS  .GT.  2)  WO  =  X(I) 

DO  140  J  =  1.  10 

120  IF  ((I.EQ.2)  .OR.  (I.EQ.4)  .OR.  (I.EQ.6))  GO  TO  130 

125  IF  (GPP(VO.C.K.S.NPL)  .GT.  0.0D0)  GO  TO  135 

WO  =  1 .02D0  *  WO 
GO  TO  125 

130  IF  (GPP(WO.C.K.N.NPL)  .LT.  0.0D0)  GO  TO  135 

WO  =  1.02D0  *  WO 
GO  TO  130 

135  STEP  =  (GREFP(WO.GS)  -  GP(WO.C.K.N.NPL))  / 

2  (GREFPP(WO.GS)  -  GPP(WO.C.K.N.NPL) ) 

IF  (STEP  .LT.  -0.03D0)  STEP  =  -0.03D0 
IF  (STEP  .GT.  +0.03D0)  STEP  =  +0.03D0 
W1  =  WO  -  STEP 

Z  =  GREFP (  Wl.GS  )  -  GP(  Wl.C.K.N.NPL  ) 

IF  (DABS(  Z  )  .LT.  0.00005D0)  GO  TO  145 
WO  =  W1 
I FLAG  =  .TRUE. 

140  CONTINUE 

145  X(I)  =  W1 

150  CONTINUE 

C  HERE  ON  TERMINATION  OF  LOOP  TO  FIND  TRUE  CRITICAL  FREQUENCIES. 

IF  (10  .GT.  0)  WRITE  (6.940)  (X(I).I  =  1,  N3) 

940  FORMAT  (  /  '  RE-COMPUTED  CRITICAL  FREQUENCIES  /  7F10.5  ) 

PASS  =  PASS  +  1 

IF  (.NOT.  I FLAG)  GO  TO  160 

IF  (PASS  .LT.  10)  GO  TO  35 
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C  JUMP  HERE  AFTER  COMPLETING  FINAL  PASS  THROUGH  SYNTHESIS. 

160  N2  =  N  +  1 

C  LOAD  DENOMINATOR  COEFFICIENTS  INTO  OUTPUT  VECTOR. 

DO  165  1=1.  N2 

XX(I)  =  C(N2  +1-1) 

165  CONTINUE 

C  LOAD  CRITICAL  FREQUENCIES  INTO  OUTPUT  VECTOR. 

DO  170  1=1,  N3 

XX(N2  +  I)  =  X(I) 

170  CONTINUE 

IF  (10  .GT.  0  )  WRITE  (6.980)  (XX(I),  1=1.  N+l  ) 

980  FORMAT  (/  *  DENOMINATOR  COEFFICIENTS:  (Constant  first)'  /  7F11.6  ) 

IF  (10  .LT.  0)  RETURN 

C  FIND  MAX  VALUE  OF  GAIN  FUNCTION,  SO  WE  CAN  NORMALIZE  TO  1.0. 

CALL  MXGAIN(  GMAX.X.C.K.N.NPL  ) 

C  NORMALIZE  GAIN  FUNCTION,  AND  COMPUTE  RESPONSE. 

WRITE  (6.950) 

950  FORMAT  (  /  T17 .  'GAIN*  /  T7.  'FREQ',  T17 .  *(DB)'  ) 

XINC  =  (X(N3)  -  X(l))  /  10.0D0 
DO  180  I  =  1.  11 
C  IN  BAND  RESPONSE. 

W1  =  X(l)  +  FLOAT(I-l)  *  XINC 

G1  =  10.0D0*DLOG10(  GAIN2(  Wl.C.K.N.NPL  )  /  GMAX  ) 

W2  =  DSQRT(  W1  ) 

WRITE  (6,960)  W2,  Gi 

960  FORMAT  (  T4.  F7.3.  F10.3  ) 

180  CONTINUE 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GAIN 2  (  W,  C.  K.  N,  NPL  ) 

DOUBLE  PRECISION  W,  C(7),  Z 
INTEGER  N,  K.  NPL 
Z  =  W  *  C(l) 

DO  10  I  =  2.  N 

Z  =  W  *  (Z  +  C(I)) 

10  CONTINUE 

Z  =  Z  +  C(N+1) 

IF  (K  .EQ.  0)  GO  TO  20 

IF  (NPL  .EQ.  0)  GAIN 2  =  (W**K)  /  Z 

IF  (NPL  .NE.  0)  GAIN2  =  (W**K)*(W+1 .D0)**NPL  /  Z 

RETURN 

20  IF  (NPL  .EQ.  0)  GAIN 2  =  1  .DO  /  Z 

IF  (NPL  .NE.  0)  GAIN2  =  ( 1 .D0+W)**NPL  /  Z 
RETURN 
END 

SUBROUTINE  MXGAIN(  GMAX,  X,  C,  K,  N,  NPL  ) 

C  FIND  MAXIMUM  GAIN  POINT  IN  AND  ABOVE  PASSBAND. 

DOUBLE  PRECISION  X(7).  C(7),  W,  GMAX,  XINC,  G 
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INTEGER  N.  K.  NPL 

XINC  =  (X(N+1)  -  X(l))  /  40.0D0 

IF  (N  .EQ.  3  .OR.  N  .EQ.  5)  XINC  =  (X(N)-X(l))  /  40.D0 
W  =  X(  1) 

GMAX  =  -25.0D0 
DO  10  I  =  1.  80 
W  =  W  +  XINC 

G  =  GAIN2(  W,  C.  K.  N.  NPL  ) 

IF  (G  .GT.  GMAX  )  GMAX  =  G 
CONTINUE 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GP  (  W,  C,  K,  N.  NPL  ) 

COMPUTE  VALUE  OF  FIRST  DERIVATIVE  OF  GAIN  FUNCTION  AT  W. 

DOUBLE  PRECISION  W.  C(7),  W1 .  W2 
INTEGER  K.  N,  NPL 
W1  =  W  +  0.0005D0 

W2  =  W  -  0.0005D0 

GP  =  GAIN2(W1,C,K.N.NPL)/0.001D0  -  GAIN2(W2.C.K.N.NPL)/0.001D0 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GPP  (  W,  C.  K.  N.  NPL  ) 

COMPUTE  VALUE  OF  SECOND  DERIVATIVE  OF  GAIN  FUNCTION  AT  W. 
DOUBLE  PRECISION  W.  C(7),  W1 ,  W2 
INTEGER  K.  N,  NPL 
W1  =  W  +  0.0005D0 

W2  =  W  -  0.0005D0 

GPP  =  GP(VI.C,K.N.NPL)/0.001D0  -  GP(W2 .C.K ,N,NPL)/0 .001D0 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GREF  (  W,  GS  ) 

COMPUTE  REFERENCE  GAIN  AT  W. 

DOUBLE  PRECISION  W,  GS 

GREF  =  DEXPC  0.38245111D0  *  GS  *  DLOGIOC  W  )  ) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  GREFP  (  W.  GS  ) 

COMPUTE  FIRST  DERIVATIVE  OF  GREF(X)  AT  W. 

DOUBLE  PRECISION  W,  GS 

GREFP  =  0. 166096 4D0  *  GS  *  GREF(  W,  GS  )  /  W 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  GREFPP  (  W,  GS  ) 

COMPUTE  SECOND  DERIVATIVE  OF  GREF(X)  AT  W. 

DOUBLE  PRECISION  W.  GS 

GREFPP  =  (GS  *  0.166096 ADO  /  W)*(  GREFP(W.GS)  -  GREF(W,GS)/W  ) 
RETURN 


to  O 


128 


DL1HITL 


A.l  DLIMTL 

DOUBLE  PRECISION  FUNCTION  DLIMIT  (  RIP,  BW.  GS.  K,  X2  ) 

DOUBLE  PRECISION  RIP.  BW,  GS.  XI.  X2,  Gl,  G2.  G3.  SG22 
XI  =  (1. DO  /  BW)  **  2 

CALL  SECANT  (  XI.  RIP.  GS.  K.  X2  ) 

Gl  =  DEXP  (0.3824511D0  *  GS  *  DLOGIO(XI)  -  RIP  *  0.46051702D0  ) 
G2  =  DEXP  (0.3824511D0  *  GS  *  DLOG10CX2) ) 

G3  =  DEXP  (  -RIP  *  0.46051702DO  ) 

SG22  =  0 .166096 4D0  *  GS  /  (G2  *  X2  ) 

T1  =  DNUMC  X1.X2.G1.G2.G3.SG22.K  ) 

T2  =  DDEN  (  X1.X2.K  ) 

DLIMIT  =  T1  /  T2 

RETURN 

END 


SUBROUTINE  SECANT  (  XI.  RIP.  GS,  K.  X2  ) 

DOUBLE  PRECISION  INC.  TEMP,  B1 .  B2.  B2N.  X2L,  X2H.  X2N. 

2  XI.  RIP.  GS.  X2 

B1  =  FUNCT  (  XI.  XI.  RIP,  GS,  X  ) 

X2L  =  XI 

INC  =  (1.D0  -  XI)  /  10 .DO 
X2H  =  X2L  +  INC 

C  Search  for  an  argument  that  produces  a  sign  change  from  that  of  B1 . 

DO  10  I  =  1,  10 

B2  =  FUNCTC  XI.  X2H.  RIP,  GS,  K  ) 

IF  ( B2  .EQ.  0  .DO)  GO  TO  70 

IF  (B2  .GT.  0.D0  .AND.  B1  .LT.  0.D0)  GO  TO  20 

IF  (B1  .GT.  0.D0  .AND.  B2  .LT.  0.D0)  GO  TO  20 

X2L  =  X2H 
B1  =  B2 

X2H  =  X2H  +  INC 
10  CONTINUE 

X2L  and  X2H  nov  bracket  the  solution.  Get  improved  estimate  of  solution. 
0  TEMP  =  -1.2D0  *  B1  /  B2 

IF  (TEMP  .LT.  1 .DO)  X2N  =  X2L  +  TEMP  *  (X2H  -  X2L)  /  2.D0 

IF  (TEMP  ,GT.  1 .DO)  X2N  =  X2H  -  (X2H  -  X2L)  /  (2.D0  *  TEMP) 

B2N  =  FUNCTC  XI.  X2N,  RIP.  GS,  K  ) 

IF  (  DABS(B2N)  .LT.  1 .OD-7  )  GO  TO  60 
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IF  (B2N  .LT.  O.DO)  GO  TO  4 0 
IF  (B1  .GT.  O.DO)  GO  TO  30 
B2  =  B2N 
X2H  =  X2N 
GO  TO  20 

30  B1  =  B2N 

X2L  =  X2N 
GO  TO  20 

40  IF  (B1  .GT.  O.DO)  GO  TO  50 

B1  =  B2N 
X2L  =  X2N 
GO  TO  20 

50  B2  =  B2N 

X2H  =  X2N 
GO  TO  20 

60  X2  =  X2N 
RETURN 

70  X2  =  X2H 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  ZERO  (  C  ) 

C  Zero  out  each  element  of  array  C(  ) . 
DOUBLE  PRECISION  C(6) 

DO  10  1=1.6 

10  C(I)  =  O.DO 

RETURN 
END 


DOUBLE  PRECISION  FUNCTION  POLY  (  C,  X  ) 

DOUBLE  PRECISION  C(6).  X 

C  C( 1)  is  constant  coefficient. 

C  5th  order  is  maximum  assumed. 

C  Value  of  polynomial  at  X  is  returned. 

POLY  =  X  *  C(6) 

DO  10  I  =  5.  2,  -1 
10  POLY  =  X  *  (  POLY  +  C(I)  ) 

POLY  =  POLY  +  C(l) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  FUNCT  (  XI,  X2.  RIP.  GS,  K  ) 

DOUBLE  PRECISION  XI .X2.RIP.GS ,C( 6) ,G1 .G2.G3 .SG22 .ALPHA, BETA, 
2  GAMMA ,  DELTA 

CALL  ZERO  (C) 


I 
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G1  =  DEXF(  0.382451 IDO  *  GS  *  DLOGIO(XI)  -  RIP  *  0.46051702D0  ) 

G2  =  DEXP(  0.3 82 451 IDO  *  GS  *  DLOG10(X2)  ) 

G3  =  DEXP(  -  RIP  *  0.46051702D0  ) 

SG22  =  0. 166096 4D0  *  GS  /  C  G2  *  X2  ) 

ALPHA  =  1.D0/G3 
BETA  =  X1**K  /  G1 
GAMMA  =  X2**K  /  G2 

DELTA  =  K*(X2**(K-1))/G2  -  SG22*(X2**K) 

C(l)  =  BETA  -  GAMMA  -  XI* (DELTA  -  DELTA*X1) 

2  ♦  (GAMMA  -  ALPHA) *(X1**2) 

C( 2)  =  2 .D0*ALPHA*X1  +  DELTA  +  2 .DO* (GAMMA  -  BETA) 

2  -  (2.D0*GAMMA*X1  +  DELTA*(X1**2)  ) 

C(3)  =  DELTA  *  XI  +  BETA  -  DELTA  -  ALPHA 

FUNCT  =  POLY  (  C,  X2  ) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  DNUM  (  XI ,X2,G1 ,G2 ,G3 .SG22.K  ) 

DOUBLE  PRECISION  XI ,X2 ,G1 ,G2 ,G3 .SG22 ,C( 6) .ALPHA. BETA .GAMMA, DELTA 
CALL  ZERO  (C) 

ALPHA  =  1.D0  /  G3 
BETA  =  X1**K  /  G1 
GAMMA  =  X2**K  /  G2 

DELTA  =  K*(X2**(K-1))/G2  -  SG22*(X2**K) 

C(l)  =  GAMMA  *  (Xl**2)  *  (XI  -  1 .DO) 

C(2)  =  Xl*(  2 .D0*GAMMA  +  DELTA*X1  -  (DELTA  +  2 .D0*GAMMA)*(X1**2) ) 
C(3)  =  -BETA  +  Xl*(  -DELTA  -  3 .D0*GAMMA  +  3 .D0*GAMMA*X1  + 

2  (DELTA  +  ALPHA)*(X1**2)  ) 

C( 4)  =  2.D0*BETA  +  DELTA*X1  -  (DELTA  +  2 .D0*ALPHA)*X1**2 
C( 5)  =  ALPHA  *  XI  -  BETA 
DNUM  =  POLY  (  C,  X2  ) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  DDEN  (  XI,  X2.  K  ) 

DOUBLE  PRECISION  C(6),X1,X2 
CALL  ZERO  (C) 

C(l)  =  (Xl**2)  *  (  XI  -  1.D0  ) 

C( 2)  =  2. DO  *  XI  *  (  1.D0  -  (Xl**2)  ) 

C(3)  =  Xl**3  +  3.D0*(X1**2)  -  3 .D0*X1  -  1 .DO 
C( 4)  =  2. DO  *  (1. DO  -  Xl**2) 

C( 5)  =  XI  -  1 .DO 
DDEN  =  POLY(  C,  X2  ) 

RETURN 

END 
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L.k  DLIMITD 

DOUBLE  PRECISION  FUNCTION  DLIMIT  C  RIP,  XI,  X3.  GS,  K,  Q.  X2  ) 

DOUBLE  PRECISION  RIP.  GS .  Xl ,  X2.  X3.  G1 ,  G2.  G3 ,  SG22 
INTEGER  Q 

CALL  SECANT  (  XI,  X3.  RIP,  GS,  K,  Q.  X2  ) 

G1  =  DEXP  (0.382451 IDO  *  GS  *  DLOGIO(XI)  -  RIP  *  0.4605I702D0  ) 

G2  =  DEXP  (0.3824511D0  *  GS  *  DLOG10(X2)) 

G3  =  DEXP  (0.382451 IDO  *  GS  *  DL0G10(X3)  -  RIP  *  0.46O51702D0  ) 

SG22  =  0 .166096 4D0  *  GS  /  (G2  *  X2  ) 

T1  =  DNUM(  X1,X2,X3,G1,G2,G3.SG22,K,Q  ) 

T2  =  DDEN  (  X1.X2.X3  ) 

DLIMIT  =  T1  /  T2 

RETURN 

END 


SUBROUTINE  SECANT  (  XI.  X3.  RIP,  GS.  K.  Q,  X2  ) 

DOUBLE  PRECISION  INC,  TEMP,  B1 .  B2.  B2N.  X2L.  X2H,  X2N, 

2  XI.  RIP.  GS.  X2.  X3 
INTEGER  Q 

B1  =  FUNCT  (  XI,  XI.  X3,  RIP,  GS,  K,  Q  ) 

X2L  =  XI 

INC  =  (X3  -  XI)  /  10. DO 
X2H  =  X2L  +  INC 

C  Search  for  an  argument  that  produces  a  sign  change  from  that  of  B1 . 

DO  10  I  =  1.  10 

B2  =  FUNCT(  XI,  X2H,  X3,  RIP,  GS ,  K,  Q  ) 

IF  ( B2  .EQ.  0 .DO)  GO  TO  70 

IF  (B2  .GT.  0 .DO  .AND.  B1  .LT.  0.D0)  GO  TO  20 

IF  ( B1  .GT.  0 .DO  .AND.  B2  .LT.  0.D0)  GO  TO  20 

X2L  =  X2H 
B1  =  B2 

X2H  =  X2H  +  INC 
10  CONTINUE 

C  X2L  and  X2H  now  bracket  the  solution.  Get  improved  estimate  of  solution. 
20  TEMP  =  -B1  /  B2 

IF  (TEMP  .LT.  1 .DO)  X2N  =  X2L  TEMP  *  (X2H  -  X2L)  /  2.D0 

IF  (TEMP  .GT.  1 .DO)  X2N  =  X2H  -  (X2H  -  X2L)  /  (2. DO  *  TEMP) 

B2N  =  FUNCT (  XI,  X2N,  X3,  RIP,  GS,  K,  Q  ) 

IF  (  DABS( B2N)  .LT.  1 .OD-7  )  GO  TO  60 
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IF  (B2N  .LT.  O.DO)  GO  TO  40 
IF  (B1  .GT.  O.DO)  GO  TO  30 
B2  =  B2N 
X2H  =  X2N 
GO  TO  20 

30  B1  =  B2N 

X2L  =  X2N 
GO  TO  20 

40  IF  (B1  .GT.  O.DO)  GO  TO  50 

B1  =  B2N 
X2L  =  X2N 
GO  TO  20 

50  B2  =  B2N 

X2H  =  X2N 
GO  TO  20 

60  X2  =  X2N 
RETURN 

70  X2  =  X2H 
RETURN 
END 

DOUBLE  PRECISION  FUNCTION  ZERO  (  C  ) 

C  Zero  out  each  element  of  array  C(  ) . 
DOUBLE  PRECISION  C(6) 

DO  10  1=1,6 

10  C(I)  =  O.DO 

RETURN 
END 


DOUBLE  PRECISION  FUNCTION  POLY  (  C,  X  ) 

DOUBLE  PRECISION  C(6),  X 

C  C( 1)  is  constant  coefficient. 

C  5th  order  is  maximum  assumed. 

C  Value  of  polynomial  at  X  is  returned. 

POLY  =  X  *  C(6) 

DO  10  I  =  5,  2,  -1 
10  POLY  =  X  *  (  POLY  +  C(I)  ) 

POLY  =  POLY  +  C(l) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  FUNCT  (  XI,  X2.  X3.  RIP,  GS,  K,  Q  ) 

DOUBLE  PRECISION  XI ,X2 ,X3 .RIP.GS ,C( 6) ,G1 ,G2,G3 ,SG22, ALPHA .BETA, 
2  GAMMA,  DELTA 

INTEGER  Q 
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CALL  ZERO  (C) 

G1  =  DEXP(  0.382 451 1D0  *  GS  *  DLOGIO(XI)  -  RIP  *  0.46051702D0  ) 

G2  =  DEXP(  0.382451 IDO  *  GS  *  DLOG10(X2)  ) 

G3  =  DEXP(  0.3824511D0  *  GS  *  DLOG10(X3)  -  RIP  *  0.46051702D0  ) 

SG22  =  0 .166096 4D0  *  GS  /  (  G2  *  X2  ) 

ALPHA  =  (  X1**K  *  ( 1 .DO  +  X1)**Q  )  /  G1 

BETA  =  (  X2**K  *  (1  .DO  +  X2)**Q  )  /  G2 

GAMMA  =  (  X3**K  *  (l.DO  +  X3)**q  )  /  G3 

DELTA  =  (K*X2**(K-1 )*( 1 .D0+X2)**Q  +  Q*X2**K*( 1 .D0+X2)**(Q-1)  )/G2 
2  -  SG22  *  (X2**K  *  (l.DO  +  X2)**Q  ) 

C(l)  =  (X1**2)*(BETA  -  GAMMA  +  X3*DELTA)  +  (X3**2)  * 

2  (ALPHA  -  BETA  -  DELTA*X1) 

C( 2)  =  2 .DO* (GAMMA  -  BETA)*X1  +  2.D0*BETA*X3  +  DELTA*(X3**2) 

2  -  2  .D0*ALPHA*X3  -  DELTA*(X1**2) 

C(3)  =  DELTA  *  (XI  -  X3)  +  ALPHA  -  GAMMA 

FUNCT  =  POLY  (  C,  X2  ) 

RETURN 

END 

DOUBLE  PRECISION  FUNCTION  DNUM  (  XI .X2.X3 ,G1 ,G2 ,G3 .SG22.K ,Q  ) 

DOUBLE  PRECISION  XI ,X2 ,G1 .G2 ,G3 .SG22 ,C( 6) .ALPHA .BETA. GAMMA. DELTA. 
2  X3 


CALL  ZERO  (C) 

ALPHA  =  (  X1**K  *  (l.DO  *  X1)**Q  )  /  G1 

BETA  =  (  X2**K  *  (l.DO  +  X2)**Q  )  /  G2 

GAMMA  =  (  X3**K  *  (l.DO  +  X3)**Q  )  /  G3 

DELTA  =  (K*X2**(K-1)*(1 .D0+X2)**Q  +  Q*X2**K*(1 .D0+X2)**(Q-1)  )/G2 
2  -  SG22  *  (X2**K  *  (l.DO  +  X2)**Q  ) 

C(l)  =  (Xl**2)  *  (X3**2)  *  BETA  *  (XI  -  X3) 

C( 2)  =  XI  *  X3  *  (  DELTA*X1*(X3**2)  -  (X1**2)*(DELTA*X3  + 

2  2.D0*BETA)  +  2. DO  *  BETA  *  (X3**2)  ) 

C(3)  =  (Xl**3)  *  (GAMMA  +  DELTA*X3)  +  3 .D0*BETA*(X1**2)*X3 
2  -  ALPHA*(X3**3)  -  X1*(X3**2)*(DELTA*X3  +  3 .D0*BETA) 

C( 4)  =  DELTA*X1*(X3**2)  +  2 ,D0*ALPHA*(X3**2) 

2  -  (X1**2)*(2.D0*GAMMA  ♦  DELTA*X3) 

C( 5)  =  -ALPHA  *  X3  +  GAMMA  *  XI 
DNUM  =  POLY  (  C,  X2  ) 

RETURN 


END 
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DOUBLE  PRECISION  FUNCTION  DDEN  (  XI,  X2.  X3  ) 

DOUBLE  PRECISION  C< 6) ,X1 .X2.X3 
CALL  ZERO  (C) 

C(l)  =  (Xl**3)  *  (X3**2)  -  (Xl**2)  *  (X3**3) 

C( 2)  =  2. DO  *  XI  *  X3*(  X3**2  -  Xl**2  ) 

C(3)  =  Xl**3  +  3 .D0*(X1**2)*X3  -  3 .D0*X1*(X3**2)  -  X3**3 
C(  4)  =  2. DO  *  (X3**2  -  Xl**2) 

C( 5)  =  XI  -  X3 
DDEN  =  POLY(  C.  X2  ) 

RETURN 

END 


Appendix  B 

Odd-Order  Leaped  Kleeeet  Betworks 
Tables  of  Element  Telees 


This  appendix  contains  tables  of  element  values  for  several  typical 
3r<*  and  5**1  order  lumped  element  matching  network  topologies.  Examples 
are  shown  for  both  152  and  1002  (octave)  bandwidth. 

The  ripple  values  listed  are  in  dB.  capacitance  and  inductance  values 
are  in  Farads  and  Henries,  respectively,  and  resistances  are  in  Ohms.  All 
values  listed  are  for  networks  normalized  to  1  Ohm  source  impedance  and  an 
upper  band  edge  of  1  Radian  per  second.  To  obtain  element  values  for  any 
arbitrary  network  source  impedance.  R*  (in  Ohms),  and  any  arbitrary  upper 
band  edge  radian  frequency  ( u *  ~  2nf*) .  first  multiply  all  resistance  and 

it  ★ 

inductance  values  by  R  and  divide  all  capacitance  values  by  R  .  Then, 
divide  both  the  capacitance  and  inductance  values  by  u*  (Radians/ sec) . 

All  networks  were  synthesized  for  0  dB  minimum  insertion  loss,  and 
the  denominator  polynomial  constant  was  chosen  to  maximize  gain-bandwidth 
product  in  all  cases. 

Attempts  were  made  to  synthesize  networks  with  gain  slopes  (GS)  of  0. 
3.  and  6  dB  per  octave  for  each  topology.  For  those  topologies  where  some 
combinations  of  ripple  and  gain  slope  were  unusable,  due  to  gain  overshoot 
or  other  problems,  only  those  cases  which  could  be  synthesized  with  0  dB 
minimum  insertion  loss  are  included  in  the  tables. 
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Odd-Order  Distributed  Element  Networks 
Tables  of  Kleaeat  Values 


This  appendix  contains  tables  of  element  values  for  several  typical 
3r<*  and  5**1  order  distributed  element  matching  network  topologies.  Exam¬ 
ples  are  shown  for  both  15Z  and  100Z  (octave)  bandwidth. 

The  ripple  values  listed  are  in  dB;  the  output  impedance  and  all 

line  impedances  are  in  Ohms*  assuming  a  normalised  1  Ohm  network  source 

•  ★  .  , 

impedance.  To  obtain  values  for  some  other  source  impedance*  R  (in 
Ohms)*  simply  multiply  the  table  entries  by  R*.  In  each  case  it  is 
assumed  that  all  lines  are  of  the  same  length*  and  that  this  length  is 
electrically  one-eighth  wavelength  at  the  upper  edge  of  the  passband. 

All  networks  were  synthesized  for  0  dB  minimum  insertion  loss*  and 
the  denominator  polynomial  constant  was  chosen  to  maximize  gain-bandwidth 
product  in  all  cases. 

Attempts  were  made  to  synthesize  networks  with  gain  slopes  (GS)  of  0* 
3*  and  6  dB  per  octave  for  each  topology.  For  those  topologies  where  some 
combinations  of  ripple  and  gain  slope  were  unusable*  due  to  gain  overshoot 
or  other  problems*  only  those  cases  which  could  be  synthesized  with  0  dB 
minimum  insertion  loss  are  included  in  the  tables. 
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Appendix  D 


Single-Freqwemcy  Matching  letwork  Design 


For  narrowband  amplifier  designs  (10Z  bandwidth  or  less)  it  is  often 
possible  to  obtain  satisfactory  performance  from  the  impedance  matching 
networks  if  they  are  simply  designed  to  achieve  the  required  match  at  a 
single  frequency.  This  technique  has  the  advantage  of  simplicity  (the 
sophisticated  computer  programs  required  for  broadband  synthesis  are  not 
needed) •  and  the  resulting  networks  require  the  fewest  possible  circuit 
elements.  This  technique  cannot  produce  a  good  match  over  a  wide  band,  of 
course,  nor  can  it  produce  gain  responses  of  arbitrary  slope;  performance 
in  both  of  these  areas  can  usually  be  improved  somewhat  by  computer 
optimization,  however. 

U*1  General  Deaign  Strategy 


The  FET  matching  problem  can  be  viewed  as  one  of  matching  a  capaci¬ 
tive  admittance  (the  input  or  output  of  an  FET)  to  a  resistive  impedance 

(either  the  source  or  load  for  the  amplifier).  We  wish  to  use  a  shunt 

inductor  as  the  first  element  connected  to  the  FET  (for  bias  insertion), 

with  the  rest  of  the  circuit  topology  dictated  by  the  desire  to  use  the 
fewest  possible  elements  to  complete  the  match. 

There  are  two  classes  of  solution  to  this  problem;  the  magnitude  of 
the  real  part  of  the  FET  admittance  controls  which  of  the  two  solutions 
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may  be  used  in  a  given  case.  The  methods  used  to  form  each  of  the  two 
types  of  solutions  will  be  illustrated  by  examples  in  the  following  sec¬ 
tions  . 


D. .1.1  £-1  Solution  Claflfl 

When  the  real  part  of  the  normalized  FET  admittance  is  less  than  1»  a 
simple  C-L  network  may  be  used  to  match  the  device  to  a  resistive  source 
or  load.  Figure  D.l  graphically  illustrates  the  solution  method  on  the 
Smith  Chart. 

The  FET  admittance  is  represented  at  A,  and  we  move  to  point  B  by 
adding  a  shunt  inductor.  The  size  of  the  inductor  is  chosen  so  that  when 
we  reflect  point  B  into  the  impedance  plane,  by  moving  from  B  to  C,  point 
C  falls  on  the  r  =  1  circle.  Then  the  match  is  completed  by  adding  a 
series  capacitor  to  move  from  C  to  D. 

Figure  D.2  shows  the  single-frequency  C-L  networks  that  result  from 
this  procedure  which  can  be  used  either  for  input  matching  (Figure  D.2(a)) 
or  output  matching  (Figure  D.2(b)). 
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When  the  reel  pert  of  the  normalized  FET  admittance  ie  greater  than 
1 >  en  L-C-L  network  cen  be  need  to  accompliah  the  match.  Figure  D.3 
graphically  illustrates  the  solution  method  on  the  Smith  Chart* 

The  FET  admittance  is  again  represented  at  point  A,  and  we  move  to 
point  B  by  adding  a  shunt  inductor.  This  inductor  must  be  sized  so  that  B 
t  liea  below  the  zero  ausceptance  axis,  but  otherwise  its  value  is  not  crit- 

f  J 

ical  (this  constraint  fixea  an  upper  limit  on  posaible  values  for  this 


Figure  D.3.  Smith  Chart  L-C-L  Network  Solution 

inductor)*  We  then  reflect  !  to  C  to  move  into  the  impedance  plane*  and 
then  a  aeriea  capacitor  ia  added  to  move  from  C  to  D.  Thia  capacitor  muat 
be  aized  ao  that  when  we  reflect  point  D  into  the  admittance  plane  by  mov~ 


ing  from  D  to  £•  point  E  falls  on  the  g  =  1  circle*  Then  the  match  is 
completed  by  adding  another  shunt  inductor  to  move  from  E  to  F. 

Figure  0.4  shows  the  single-frequency  L-C-L  networks  that  result  from 
this  procedure  which  can  be  used  either  for  input  siatching  (Figure  D.4(a) 
or  output  matching  (Figure  0.4(b)) . 


(a) 


Figure  D.4.  L-C-L  Input/Output  Networks 
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U..Z  Inter a tag c  networks 

Interstage  networks  can  also  be  designed  by  a  similar  Smith  Chart 
technique.  First*  a  shunt  inductor  value  is  determined  that  tranaforms 
the  input  admittance  of  the  second  FET  to  a  pure  conductance  (or  resis¬ 
tance  in  the  impedance  plane) .  Then  the  output  matching  problem  for  the 
first  FET  is  normalized  to  this  level*  and  a  C-L  or  L-C-L  network*  as 
required*  is  designed  to  complete  the  match  to  the  resistive  value  now 
presented  at  the  input  to  the  second  FET. 

The  two  circuits  are  then  combined  into  a  single  L-C-L  network  to 
form  the  interstage  matching  network.  (If  the  FET  output  problem  required 
an  L-C-L  network*  the  output  inductor  of  this  network  is  combined  with  the 
shunt  inductor  used  initially  at  the  input  to  the  second  FET  to  form  a 
single  inductor  at  that  point  in  the  circuit.) 

Figure  D.5  shows  the  L-C-L  interstage  network  topology  that  always 
results  from  this  procedure. 


Figure  D.5.  Single-Frequency  Interstage  Network 
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tt.l  Conversion  to.  Distributed  inductor* 

In  tome  circumstance*  it  may  be  desirable  to  use  shunt  short- 
circuited  transmission  line  stubs  in  lieu  of  the  lumped  inductors  used  in 
the  above  single~frequency  solutions.  A  simple  procedure  may  be  used  to 
make  this  conversion,  usually  computed  at  band  center.  At  the  chosen  fre¬ 
quency  (f).  compute  the  inductive  reactance  (Z^)  of  the  lumped  inductor 

(L)  to  be  replaced 


ZL  =  2«fl 

then  choose  a  line  length  (6.  in  wavelengths)  and  characteristic  impedance 
(Zq)  for  the  shorted  stub  such  that 

ZL  =  Z0  tan  e- 

Over  narrow  bands  the  performance  of  the  circuits  after  conversion  to 
shunt  shorted  stubs  will  be  virtually  the  same  as  with  the  lumped  induc¬ 
tors;  over  wider  bands,  computer  optimization  of  the  circuit  may  be 
necessary. 

H.A  Computer  Implementation 

Although  single-frequency  matching  problems  may  be  solved  graphi¬ 
cally.  they  may  also  be  solved  more  conveniently  by  computer.  The  S MATCH 
(Smith  Chart  MATCH)  program  has  been  developed  to  numerically  solve  the 
same  classes  of  problems  solved  graphically  above. 
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The  FORTRAN  source  code  listing  for  the  S MATCH  progrem  (end  its  sub¬ 
routines)  follows  this  section*  end  Figure  D.6  shows  example  runs  of  the 
program  for  each  of  the  three  cases:  a  C-L  network*  an  L-C-L  network*  and 
an  interstage  network.  The  program  is  interactive*  as  can  be  seen  from 
the  examples*  and  the  solutions  are  computed  and  output  almost  immedi¬ 
ately. 
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SMITH  CHART  MATCHING  PROGRAM  VI .0  JUNE  82 

PROBLEM  TYPE  ( 0 =INPUT/OUTPUT .  1 ^INTERSTAGE )  : 

0 

FET  MODEL  VALUES:  R(OHMS).  C(PF)  .  FREQ(GHZ)  . 

MODEL  TYPE  (1=SERIES.  0=SHUNT) 

9.268  0.5099  11  .6  1 

MATCHING  NETWORK  AT  11.60  GHZ  : 

50  OHM  SOURCE 

SERIES  CAPACITOR  :  .317  PF 

SHUNT  INDUCTOR  :  .318  NH 

FET  :  R  =  9.268  (SERIES  MODEL) 

C  =  .510 

PROBLEM  TYPE  ( 0 =INPUT/OUTPUT  .  1 =INTERSTAGE )  : 

0 

FET  MODEL  VALUES:  R(OHMS).  C(PF),  FREQ(GHZ). 

MODEL  TYPE  (1=SERIES.  0=SHUNT) 

16  .32  0  .3654  11  .6  0 

L-C-L  NETWORK  REQUIRED.  CHOOSE  L2 .  .52  MAX  ALLOWED. 

0  .396 

MATCHING  NETWORK  AT  11.60  GHZ  : 

50  OHM  SOURCE 

SHUNT  INDUCTOR  :  .472  NH 

SERIES  CAPACITOR  :  .539  PF 

SHUNT  INDUCTOR  :  .396  NH 

FET  :  R  =  16.320  (SHUNT  MODEL) 

C  =  .365 

PROBLEM  TYPE  ( 0 =INPUT/OUTPUT  ,  1 = INTERSTAGE )  : 

1 

MODEL  VALUES  FOR  OUTPUT  OF  FIRST  FET  : 

R(OHMS).  C(PF),  FREQ(GHZ),  MODEL  TYPE  (1=SERIES.  0=SHUNT) 
199.5  0.1058  11.6  0 

MODEL  VALUES  FOR  INPUT  TO  SECOND  FET  : 

R(OHMS).  C ( PF ) •  MODEL  TYPE  (1=SERIES.  0=SHUNT) 

9.268  0.5099  1 

INTERSTAGE  MATCHING  NETWORK  AT  11.60  GHZ  : 

OUTPUT  OF  FIRST  FET  :  R  =  199.500  (SHUNT  MODEL) 

C  =  .106 

SHUNT  INDUCTOR  :  1.025  NH 

SERIES  CAPACITOR  :  .139  PF 

SHUNT  INDUCTOR  :  .413  NH 

INPUT  OF  SECOND  FET  :  R  =  9.268  (SERIES  MODEL) 

C  =  .510 


Figure  D.6.  S MATCH  Program  Example  Rune 
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C  SMITH  CHART  MATCHING  PROGRAM  J.T.  DIJAK  21  JUNE  82 
C  FILENAME:  S MATCH  FORTRAN 

REAL  PI .FR.FC .FREQ ,FZC .FRN.FZCN ,C .LI .L2 ,L3 ,YL *FR2 ,FC2 
COMPLEX  YF.  ZF 

,  INTEGER  TYPE1 •  TYPE2 

PI  =  3.1415926 
WRITE  (6.890) 

890  FORMAT  (//•  SMITH  CHART  MATCHING  PROGRAM' ,5X. 'VI .0 ' ,5X. 'JUNE  82') 

« 

C  CHOOSE  PROBLEM  TYPE. 

5  WRITE  (6.900) 

900  FORMAT  (  /  •  PROBLEM  TYPE  (0=INPUT/OUTPUT.  1 = INTERSTAGE)  ) 

READ  (5.*.END=1000)  TYPE1 
IF  (TYPE1  .EQ.  1)  GO  TO  40 

C  SOLVE  THE  INPUT/OUTPUT  NETWORK  PROBLEM. 

WRITE  (6.910) 

910  FORMAT  (  / '  FET  MODEL  VALUES:  R(OHMS).  C(PF),  FREQ(GHZ) .  '. 

2  'MODEL  TYPE  (1=SERIES,  0=SHUNT)  •  ) 

READ  (5.*)  FR.  FC.  FREQ.  TYPE1 

C  COMPUTE  CAPACITIVE  REACTANCE  OF  FET  MODEL. 

FZC  =  500.0  /  (PI  *  FREQ  *  FC) 

C  NORMALIZE  TO  50  OHMS. 

FRN  =  FR  /  50.0 
FZCN  =  FZC  /  50.0 

IF  (TYPE1  .EQ.  1)  GO  TO  10 
C  SET  UP  FET  SHUNT  MODEL. 

YF  =  CMPLX(  1.0/FRN  .  1 .0/FZCN  ) 

GO  TO  20 

C  SET  UP  FET  SERIES  MODEL. 

10  ZF  =  CMPLX(  FRN,  -FZCN  ) 

YF  =  1.0  /  ZF 

C  DETERMINE  IF  C-L  OR  L-C-L  NETWORK  IS  REQUIRED. 

20  IF  (REAL(YF)  .GT.  1.0)  GO  TO  30 

*  C  SOLVE  C-L  NETWORK  FOR  ELEMENT  VALUES. 

CALL  CLNET  (  YF.  FREQ,  50.0,  C,  LI  ) 

C  PRINT  SOLUTION. 

WRITE  (6.920)  FREQ.  C,  LI 

„  920  FORMAT  (  /  '  MATCHING  NETWORK  AT  ',  F5.2,  '  GHZ  :'  / 

2  T5.  '50  OHM  SOURCE'  / 

3  T7 ,  'SERIES  CAPACITOR  :  ',  F7.3.  '  PF'  / 

4  T7 ,  'SHUNT  INDUCTOR  :  ',  F7.3.  '  NH'  ) 

IF  (TYPE1  .EQ.  1)  WRITE  (6.980)  FR.  FC 

IF  (TYPE1  .EQ.  0)  WRITE  (6.982)  FR,  FC 

980  FORMAT (  T5.  'FET  :  R  =  ' »  F8.3,  5X,  '(SERIES  MODEL)'  / 
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2  T13»  ’C  =  ’.  F8.3  ) 

982  FORMAT  (  T5.  ’FET  :  R  =  ’ »  F8.3.  5X.  ’(SHUNT  MODEL)’  / 

2  T13.  ’C  =  ’.  F8.3  ) 

GO  TO  5 

C  SOLVE  L-C-L  NETWORK  FOR  ELEMENT  VALUES. 

30  CALL  LCLNET(  YF.  FREQ.  50.0.  LI.  C.  L2  ) 

C  PRINT  SOLUTION. 

WRITE  (6.940)  FREQ.  LI.  C.  L2 
940  FORMAT  (  /’  MATCHING  NETWORK  AT  ’.  F5.2,  ’GHZ  :'  / 

2  T5.  *50  OHM  SOURCE’  / 

3  T7 .  ’SHUNT  INDUCTOR  :  ’.  F7.3.  '  NH’  / 

4  T7 .  ’SERIES  CAPACITOR  :  ’.  F7.3.  ’  PF’  / 

5  T7 .  ’SHUNT  INDUCTOR  :  ’.  F7.3.  '  NH’  ) 

IF  (TYPEl  .EQ.  1)  WRITE  (6.980)  FR.  FC 

IF  (TYPEl  .EQ.  0)  WRITE  (6.982)  FR.  FC 

GO  TO  5 

C  SOLVE  INTERSTAGE  PROBLEM. 

40  WRITE  (6.950) 

950  FORMAT  (/'  MODEL  VALUES  FOR  OUTPUT  OF  FIRST  FET  ;•  / 

2  5X. ’R(OHMS) .  C(PF) •  FREQ(GHZ) .  MODEL  TYPE  (1=SERIES.  0=SHUNT)’) 

READ  (5.*)  FR.  FC,  FREQ.  TYPEl 
WRITE  (6.960) 

960  FORMAT  (/'  MODEL  VALUES  FOR  INPUT  TO  SECOND  FET  : •  / 

2  5X,  'R(OHMS) .  C(PF).  MODEL  TYPE  (1=SERIES.  0=SHUNT)  •  ) 

READ  (5.*)  FR2.  FC2,  TYPE2 

C  MODEL  THE  INPUT  TO  THE  SECOND  FET  AS  AN  ADMITTANCE  YF. 

C  COMPUTE  CAPACITIVE  REACTANCE  OF  MODEL. 

FZC  =  500.0  /  (PI  *  FREQ  *  FC2) 

IF  (TYPE2  .EQ.l)  GO  TO  50 
C  SHUNT  MODEL  FOR  FET  2. 

YF  =  CMPLX(  1.0/FR2.  1 .0/FZC  ) 

GO  TO  60 

C  SERIES  MODEL  FOR  FET  2. 

50  ZF  =  CMPLX(  FR2.-FZC  ) 

YF  =  1.0  /  ZF 

C  RESONATE  OUT  SUSCEPTANCE  OF  YF  USING  SHUNT  INDUCTOR  L2. 

60  YL  =  AIMAG(  YF  ) 

L2  =  0.5  /  (PI  *  FREQ  *  YL) 

C  COMPUTE  NORMALIZATION  CONSTANT  FOR  THE  REST  OF  THE  PROBLEM. 

RN  =  1.0  /  REAL(  YF  ) 

C  SET  UP  NORMALIZED  MODEL  FOR  OUTPUT  OF  FIRST  FET. 

FZC  =  500.0  /  (PI  *  FREQ  *  FC) 
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C  NORMALIZE  MODEL  VALUES. 

FRN  =  FR  /  RN 
FZCN  =  FZC  /  RN 

IF  (TYPE1  .EQ.  1)  GO  TO  70 
C  SHUNT  MODEL. 

YF  =  CMPLX(  I.O/FRN,  1 .O/FZCN  ) 

GO  TO  80 

C  SERIES  MODEL . 

70  ZF  =  CMPLX(  FRN,  -FZCN  ) 

YF  =  1.0  /  ZF 

C  DETERMINE  IF  C-L  OR  L-C-L  NETWORK  IS  REQUIRED. 

80  IF  (REAL(YF)  .GT.  1.0)  GO  TO  100 

C  SOLVE  C-L  NETWORK  FOR  ELEMENT  VALUES. 

CALL  CL NET (  YF.  FREQ,  RN.  C,  LI  ) 

C  PRINT  SOLUTION. 

90  WRITE  (6.970)  FREQ 

970  FORMAT  (  /  1  INTERSTAGE  MATCHING  NETWORK  AT  ’,  F5.2.  * 
IF  (TYPE1  .EQ.  1)  WRITE  (6,972)  FR.  FC 

IF  (TYPEl  .EQ.  0)  WRITE  (6,974)  FR.  FC 

972  FORMAT  (  T5.  ’OUTPUT  OF  FIRST  FET  :  R  =  ',  F8.3.  5X. 

2  ’(SERIES  MODEL)’  /  T29,  *C  =  '.  F8.3  ) 

974  FORMAT  (  T5.  ’OUTPUT  OF  FIRST  FET  :  R  =  '.  F8.3,  5X. 
2  ’(SHUNT  MODEL)’  /  T29.  ’C  =  ',  F8.3  ) 

WRITE  (6.976)  LI.  C,  L2 

976  FORMAT  (  T7 .  ’SHUNT  INDUCTOR  :  ’.  F7.3,  ’  NH'  / 

2  T7 ,  ’SERIES  CAPACITOR  :  F7.3,  '  PF*  / 

3  T7 ,  ’SHUNT  INDUCTOR  :  F7.3,  »  NH’  ) 

IF  (TYPE2  .EQ.  1)  WRITE  (6,978)  FR2 ,  FC2 

IF  (TYPE2  .EQ.  0)  WRITE  (6.979)  FR2,  FC2 

978  FORMAT  (  T5.  'INPUT  OF  SECOND  FET  :  R  =  ’,  F8.3.  5X, 

2  '(SERIES  MODEL)’  /  T29.  ’C  =  ’,  F8.3  ) 

979  FORMAT  (  T5.  ’INPUT  OF  SECOND  FET  :  R  =  ’,  F8.3*  5X. 

2  '(SHUNT  MODEL)'  /  T29,  'C  =  *,  F8.3  ) 

GO  TO  5 

C  SOLVE  L-C-L  NETWORK  PROBLEM. 

100  CALL  LCLNET(  YF,  FREQ.  RN.  L3.  C,  LI  ) 

C  FORM  PARALLEL  COMBINATION  OF  L2  &  L3 . 

L2  =  L2  *  L3  /  (L2  +  L3) 

C  PRINT  SOLUTION . 

GO  TO  90 


GHZ 


1000  STOP 
END 


o  o 
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SUBROUTINE  SOLVE  (  YF,  YL  ) 

C  FIND  AN  INDUCTIVE  SUSCEPTANCE  VALUE.  YL,  SUCH  THAT  THE 
C  NORMALIZED  SMITH  CHART  MATCHING  PROBLEM  IS  SOLVED  - 
C  FROM  A  CAPACITIVE  ADMITTANCE  STARTING  POINT.  YF. 

C  SOLUTION  IS  GUARANTEED  FOR  ALL  CAPACITIVE  YF. 

C  THE  PROBLEM  IS  SOLVED  BY  BISECTION,  WITH  THE  FINAL 

C  RELATIVE  ERROR  IN  YL  LESS  THAN  0.1Z  . 

C  ON  ENTRY.  YF  IS  THE  CAPACITIVE  ADMITTANCE  STARTING  POINT. 

C  ON  EXIT  .  YL  IS  THE  SOLUTION  SUSCEPTANCE. 

C  J.T.  DIJAK  JUNE  1982  //  PART  OF  SMATCH  PROGRAM. 

REAL  YL1 ,  YL2,  YL,  F,  X,  Y,  T 
COMPLEX  YF 

C  GENERATE  AN  INITIAL  GUESS  FOR  YL  BASED  UPON  REAL(YF) .  Y  IS  THE 
C  IMAGINARY  PART  OF  THE  ADMITTANCE  OF  THE  INITIAL  GUESS. 

X  =  YF 

IF  (X  .LT.  0.4)  Y  =  0.2  +  0.3  *  X  /  0.4 
IF  (X  .GE.  0.4  .AND.  X  .LE.  0.5)  Y  =  0.5 

IF  (X  .GT.  0.5  .AND.  X  .LE.  0.95)  Y  =  0.3  +  0.3  *  (0.95  -  X)/0.45 

IF  (X  .GT.  0.95)  Y  =  (1.0  -  X)  *  0.25  /  0.05 
C  GENERATE  INITIAL  GUESS  INDUCTIVE  SUSCEPTANCE. 

YL1  =  Y  +  AIMAG(  YF  ) 

C  WORST  CASE  ERROR  IN  THIS  INITIAL  GUESS  IS  +/-  0.1. 

C  CHOOSE  A  SECOND  STARTING  POINT  FOR  OPPOSITE  SIGN  OF  ERROR  FUNCTION. 
F  =  1.0  -  REAL(  1.0  /  (YF  -  CMPLX(0 .0 ,YL1 )  )  ) 

YL2  =  YL1 

IF  (F  .GT.  0.0)  YL2  =  YL2  -  0.15 

IF  (F  .LE.  0.0)  YL2  =  YL2  +  0.15 

SOLUTION  INTERVAL  IS  NOW  DEFINED.  INSURE  THAT  YL1  CORRESPONDS 
TO  F  >0. 

IF  (F  .GT.  0.0)  GO  TO  10 
T  =  YL1 
YL1  =  YL2 

YL2  =  T 

C  BISECT  THE  INTERVAL. 

*  10  YL  =  (  YL1  +  YL2  )  /  2.0 

C  EVALUATE  THE  ERROR  FUNCTION. 

F  =  1.0  -  REAL(  1.0  /  (  YF  -  CMPLX(  0.0,  YL  )  )  ) 

%  C  STOP  ONCE  THE  TERMINATION  CRITERIA  IS  SATISFIED. 

IF  (  (YL1  -  YL2)  .LT.  (YL/1000)  )  RETURN 

C  REPLACE  ONE  OF  THE  SOLUTION  INTERVAL  BOUNDARIES  &  ITERATE. 

IF  (  F  .GT.  0.0  )  YL1  =  YL 
IF  (  F  .LT.  0.0  )  YL2  =  YL 
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GO  TO  10 
END 

SUBROUTINE  CLNET  (  YF.  FREQ,  RN,  C,  L  ) 

SOLVE  THE  C-L  NETWORK  PROBLEM.  WHERE  THE  NETWORK  IS  : 

SERIES  C,  SHUNT  L,  FET. 

ON  ENTRY,  YF  IS  STARTING  POINT  ADMITTANCE,  FREQ  IS  IN  GHZ,  AND 
RN  IS  NORMALIZATION  CONSTANT. 

ON  EXIT,  C  AND  L  ARE  THE  NETWORK  ELEMENT  VALUES  IN  PF  AND  NH. 

REAL  C,  L,  YL.  ZC,  FREQ.  RN 
COMPLEX  YF 

PI  =  3.1415926 

SOLVE  FOR  REQUIRED  SHUNT  INDUCTIVE  ADMITTANCE.  YL . 

CALL  SOLVE  (  YF,  YL  ) 

COMPUTE  SIZE  OF  SHUNT  INDUCTOR  (NH)  . 

L  =  RN  /  (2.0  *  PI  *  FREQ  *  YL) 

COMPUTE  SIZE  OF  SERIES  CAPACITOR  (PF). 

ZC  =  RN  *  AIMAG(  1.0/  (YF  -  CMPLX(  0.0.  YL  )  )  ) 

C  =  500.0  /  (PI  *  FREQ  *  ZC) 

RETURN 
END 

SUBROUTINE  LCLNET  (  YF.  FREQ,  RN,  LI.  C.  L2  ) 

SOLVE  THE  L-C-L  NETWORK  PROBLEM,  WHERE  THE  NETWORK  IS  : 

SHUNT  LI,  SERIES  C.  SHUNT  L2,  FET. 

ON  ENTRY.  YF  IS  STARTING  POINT  ADMITTANCE,  FREQ  IS  IN  GHZ.  AND 
RN  IS  NORMALIZATION  CONSTANT. 

ON  EXIT,  LI,  C.  AND  L2  ARE  THE  NETWORK  ELEMENT  VALUES  IN  PF  &  NH 

REAL  FREQ,  RN,  LI,  C,  L2,  ZC,  L2MAX,  YL1 ,  PI 
COMPLEX  YF,  Yl.  Z1 

PI  =  3.1415926 

C  FIND  MAX  ALLOWABLE  VALUE  FOR  THE  SHUNT  INDUCTOR  L2  (NEXT  TO  FET) 
C  THIS  VALUE  FOR  L2  WILL  MOVE  YF  JUST  TO  THE  REAL  AXIS. 

L2MAX  =  RN  /  (2.0  *  PI  *  FREQ  *  AIMAG(YF)  ) 

WRITE  (6,900)  L2MAX 

900  FORMAT  (  /  *  L-C-L  NETWORK  REQUIRED.  CHOOSE  L2.  ',  F7.2, 

2  '  MAX  ALLOWED.1  ) 

READ  (5,*)  L2 

IF  (L2  .GT.  L2MAX)  L2  =  0.95  *  L2MAX 

C  COMPUTE  RESULTING  ADMITTANCE  AFTER  ADDING  SHUNT  L2. 

Yl  =  YF  -  CMPLX(  0.0,  (RN/(2.0  *  PI  *  FREQ  *  L2)  )  ) 

Z1  =  1.0  /  Yl 

C  SOLVE  THE  MATCHING  PROBLEM  IN  THE  Z -PLANE  RATHER  THAN  Y-PLANE. 
CALL  SOLVE(  Zl,  ZC  ) 


1 


167 


C  COMPUTE  SIZE  OF  SERIES  CAPACITOR  (PF). 

C  =  500.0  /  (PI  *  FREQ  *  ZC  *  RN) 

C  COMPUTE  SIZE  OF  SECOND  SHUNT  INDUCTOR.  LI  (NH) . 
YL1  =  AIMAG(  1.0  /  (Z1  -  CMPLX(  0.0. ZC  )  )  ) 
LI  =  RN  /  (2.0  *  PI  *  FREQ  *  YL1) 


RETURN 

END 


